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ABSTRACT 


Information  warfere  is  a  cornerstone  of  Joint  Vision  2010  which  addresses  the 
future  strategic  environment  for  the  United  States.  An  integral  component  of  information 
warfare  is  the  continuing  development  of  joint  space  doctrine.  The  Joint  Warfare  System 
(JWARS)  is  a  large  scale,  systemic  simulation  being  developed  by  the  Joint  Warfare 
Systems  Office  to  aid  in  the  evaluation  of  future  joint  doctrine  and  force  structure.  The 
purpose  of  this  thesis  is  to  develop  and  demonstrate  Simulation  of  Satellite 
Communications  (SIMSATCOM),  a  high  resolution,  stochastic  simulation  of  satellite 
communications  for  evaluating  the  effectiveness  of  message  transmission  and  receipt  by 
specified  senders  and  receivers.  SIMSATCOM  is  designed  to  operate  as  a  stand-alone 
simulation,  but  may  be  adopted  as  a  high  resolution  module  for  a  large  scale  simulation 
such  as  JWARS.  The  thesis  describes  SIMSATCOM  in  detail  and  provides  analyses  of 
simulation  runs  for  different  jamming  levels  and  channel  capacities. 
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EXECUTIVE  SUMMARY 


Information  warfare  is  a  cornerstone  of  Joint  Vision  2010  which  addresses  the 
future  strategic  environment  for  the  United  States.  An  integral  component  of  information 
warfare  is  the  continuing  development  of  joint  space  doctrine.  The  Joint  Warfare  System 
(JWARS)  is  a  large  scale,  systemic  simulation  being  developed  by  the  Joint  Warfare 
Systems  OflSce  to  aid  in  the  evaluation  of  future  joint  doctrine  and  force  structure. 

The  purpose  of  this  thesis  is  to  develop  and  demonstrate  Simulation  of  Satellite 
Communications  (SIMSATCOM),  a  high  resolution,  stochastic  simulation  of  satellite 
communications  for  evaluating  the  effectiveness  of  message  transmission  and  receipt  by 
specified  senders  and  receivers.  SIMSATCOM  is  designed  to  operate  as  a  stand-alone 
simulation,  but  may  be  adopted  as  a  high  resolution  module  for  a  large  scale  simulation 
such  as  JWARS. 

SIMSATCOM  is  written  in  Visual  Basic  for  Apphcations  using  Microsoft  Excel. 
The  only  user  interface  is  prior  to  model  execution,  where  input  parameters  are  set.  After 
initialization,  the  model  selects  the  latitude  and  longitude  values  for  both  sender  and 
receiver.  Successful  communication  of  each  message  from  the  sender  to  the  receiver  is 
based  upon  1)  message  length,  which  is  a  function  of  the  type  of  message  transmitting;  2) 
whether  the  sender  and/or  receiver  groimd  station  is  being  totally  jammed,  partially 
jammed,  or  is  not  being  jammed;  and  3)  satellite  channel  availability.  There  are  ten 
stochastic  variables  in  SIMSATCOM  whose  values  are  dynamically  determined  by  Monte 
Carlo  procedures  using  appropriate  probability  distributions. 
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After  the  simulation  was  developed,  a  demonstration  case  was  run  to  show 
ejqjlicitly  how  the  simulation  works.  Sample  inputs  were  used,  and  four  messages  were 
selected  to  follow  through  one  replication  of  the  simulation.  The  case  was  run  until 
message  completion  occurred  for  each  of  the  four  messages.  Message  completion 
indicates  that  1)  a  message  was  transmitted  successfiilly  and  an  acknowledgment  message 
was  also  transmitted  successfully;  2)  a  message  was  transmitted  successfully,  but  the 
acknowledgment  message  was  not  transmitted  successfully;  and  after  a  period  of  time,  the 
sender  terminated  the  message  process  (this  maximum  time  to  transmit  a  message  is  a 
parameter  input  by  the  user  prior  to  starting  the  model);  and  3)  a  message  was  never 
transmitted  successfully  because  it  was  jammed  or  there  were  no  available  channels,  and 
the  sender  terminated  the  message  process. 

An  example  of  the  use  of  SIMSATCOM  for  analysis  purposes  is  presented.  A  2  x 
2  fectorial  e5q)eriment,  applyii^  the  analysis  of  variance,  was  run  using  channel  capacity 
and  anti-jam  capability  as  the  two  factors.  Five  replications  of  the  simulation  were  run  for 
each  case. 

The  model  is  a  baseline  high  resolution  communications  simulation  that  can  be 
e^anded  for  future  use.  The  model  could  be  modified  to  make  it  more  dynamic  in  terms 
of  satellite  locations;  instrumental  for  a  specified  area  of  interest;  making  a  more 
sophisticated  jamming  process;  or  as  a  baseline  for  other  types  of  satellite  models,  such  as 
near  earth  satellites  or  semi-synchronous  satellites. 
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1. 


INTRODUCTION 


Information  warfare  is  a  cornerstone  of  Joint  Vision  2010  which  addresses  the 
future  strategic  environment  for  the  United  States.  An  integral  component  of  information 
warfere  is  the  continuing  development  of  joint  space  doctrine.  The  Joint  Warfare  System 
(JWARS)  is  a  large  scale,  systemic  simulation  being  developed  by  the  Joint  Warfare 
Systems  Office  to  aid  in  the  evaluation  of  future  joint  doctrine  and  force  structure. 

The  motivation  for  this  thesis  came  from  discussions  with  the  JWARS  office  in 
Arlington,  Virginia.  They  were  interested  in  determining  whether  it  would  be  adequate  to 
model  satellites  aggregately  in  JWARS.  There  are  many  types  of  satellites  that  JWARS 
needs  to  consider;  however,  the  focus  of  this  thesis  is  to  build  a  communications 
simulation  that  can  be  expanded  to  answer  the  geosynchronous  satellite  portion  of  this 
concern.  (The  Milstar  satellites,  used  as  a  model  for  this  thesis,  are  in  inclined 
geostationary  orbit  approximately  22,400  nautical  miles  above  the  earth.  [Ref  17]) 

The  purpose  of  this  thesis  is  to  develop  and  demonstrate  Simulation  of  Satellite 
Communications  (SIMSATCOM),  a  high  resolution,  stochastic  simulation  of  satellite 
communications  for  evaluating  the  effectiveness  of  message  transmission  and  receipt  by 
specified  senders  and  receivers.  SIMSATCOM  is  designed  to  operate  as  a  stand-alone 
simulation,  but  may  be  adopted  as  a  high  resolution  module  for  a  large  scale  simulation 
such  as  JWARS. 

SIMSATCOM  has  many  parameters  that  can  be  modified  to  answer  several 
Hifffirent  questions.  By  adjusting  one  parameter  and  keeping  the  other  parameters 
constant,  n  x  m  fector  experiments  can  be  conducted. 


1 


Possible  questions  that  are  addressed  in  Chapter  V  are: 

“What  happens  if  the  satellite  capacity  is  increased?” 

“What  happens  if  the  anti-jam  capability  for  a  satellite  is  increased?” 

“Would  it  be  better  to  buy  more  anti-jam  capability  or  channel  capacity  for 
a  satellite  system?” 

To  answer  these  questions,  two  parameters  are  examined:  the  anti-jam  capability 
parameter  and  the  channel  capacity  parameter.  These  parameters  are  adjusted  to 
determine  if  the  number  of  good  messages  transmitted  increases  significantly  for  variations 
in  the  two  fectors.  In  Chapter  V,  a  sparse  2x2  factorial  experiment  is  presented  to 
demonstrate  how  the  model  can  be  used  to  conduct  analysis. 

A  background  of  existing  simulations  that  utilize  satellites  is  outlined  in  Chapter  II. 
A  description  of  how  SIMSATCOM  was  developed  and  how  it  functions  is  given  in 
Chapter  III.  Four  sanple  messages  are  tracked  through  the  model  to  show  how  each  of 
the  features  of  the  model  are  used  as  depicted  in  Chapter  IV.  A  small  experiment  was 
conducted  to  show  how  the  model  is  used  for  analysis  is  illustrated  in  Chapter  V.  Finally, 
recommendations  regarding  how  the  model  can  be  expanded  and  utili2ed  for  foture  use  is 
described  in  Chapter  VI. 
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11.  BACKGROUND 


In  1995  the  Secretary  of  Defense  tasked  DoD  to  conduct  a  study  to  analyze  the 
capabilities  of  current  models  and  simulations.  From  this  study  the  Deputy  Secretary  of 
Defense’s  Joint  Analytic  Model  Improvement  Program  (JAMIP)  was  established.  DoD 
developed  JAMIP  to  help  determine  how  to  more  effectively  conduct  joint  theater-level 
analysis.  Among  other  elements  of  JAMIP,  the  Joint  Warfare  System  (JWARS)  was 
conceptualized  to  satisfy  this  requirement.  [Ref  1 :  p.  1] 

JWARS  will  be  the  primary  model  used  by  the  services  for  analysis.  It  will 
effectively  model  ground,  air,  naval,  and  space  operations;  command,  control, 
communications,  computers,  intelligence,  surveillance  and  reconnaissance  (C4ISR); 
strategic  mobility;  logistics;  weapons  of  mass  destruction;  special  operations;  military 
operations  other  than  war  (MOOTW);  and  will  appropriately  represent  the  physical 
environment  and  its  effect  on  simulated  activities.  [Ref  1:  p2, 5]  It  will  operate  at  the 
theater,  campaign  and  mission  or  small-unit  level  for  joint  analysis  and  will  be  able  to 
model  strategic,  operational  and  tactical  levels  of  war  focusing  at  the  operational  level 
[Ref  1 :  p9]  However,  JWARS  will  not  be  interactive,  will  not  support  real-time  mission 
execution  or  have  any  real-world  direct  links  to  any  C4ISR  systems.  [Ref  2:  p2]  The 
current  plan  is  for  JWARS  Blocks  1, 11,  and  in  versions  to  be  developed  by 
December  1998,  December,  1999  and  December,  2000  respectively,  [Ref  1:  p9] 

In  order  to  correctly  model  space  requirements  there  is  a  need  to  determine  to 
what  extent  satellites  should  be  modeled  in  the  aggregate  for  the  theater  level  models. 
Previous  research  conducted  by  Captain  Robert  Payne  for  his  thesis  in  March  1995  reveals 
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that  many  of  the  primary  models  used  by  DoD  for  analysis  do  not  include  space  assets 
[Ref  3].  This  chapter  identifies  theater-level  models  and  some  stand  alone  space  models 
that  exist  today,  and  their  characteristics. 

A.  OTHER  MODELS 

1.  Integrated  Theater  Engagement  Model  (ITEM) 

ITEM  is  a  theater  level,  joint  forces  model  designed  to  support  joint  warfere 
analysis,  joint  mission  area  assessments,  major  regional  contingency  scenario  applications 
in  joint  wargames,  and  cost  and  operational  effectiveness  analysis,  among  others.  The 
proponent  for  ITEM  is  the  Navy  (N81/DNA/SPAWAR  3 1).  It  currently  models  land,  air, 
sea  and  undersea  warfare.  It  models  single  and  multiple  engagements,  and  many-on-many, 
force-on-force  and  theater  level  warfare.  It  is  a  two-sided,  deterministic,  Monte-Carlo 
simulation  that  uses  time  step  processing.  Human  participation  is  optional  for  decisions  in 
this  model.  [Ref  4] 

2.  Joint  Theater  Level  Simulation  (JTLS) 

JTLS  is  a  theater  level  model  designed  for  analysis,  development,  and  evaluation  of 
operations  plans  [Ref  5].  It  also  is  used  as  a  tool  for  force  mix  and  resource  analysis  and 
as  a  high  resolution  wargame.  JTLS  models  land,  air,  limited  naval  operations, 
intelligence  and  logistics.  It  is  a  multi-sided  model  that  requires  human  participation  for 
decisions  and  processes.  JTLS  can  model  from  a  single  brigade  to  theater  level  warfiire.  It 
is  a  dynamic,  event-step  model  that  represents  groxmd  combat  attrition  determmistically, 
and  air  and  naval  attrition  stochastically.  JTLS  requires  human  participation  for  decisions 
and  processes.  The  proponent  for  JTLS  is  the  Joint  Warfere  Center.  [Ref  6] 
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3.  Mission  Eifectiveness  Model  (MEM) 

MEM  is  a  two-sided  ballistic  missile  defense  and  ground-based  Anti-Satellite 
(ASAT)  analysis  model.  The  proponent  for  MEM  is  the  Headquarters,  United  States  Air 
Force,  Space  Division.  It  focuses  on  nuclear  ballistic  missiles  for  offensive  forces  and 
non-nuclear  weapons  (e.g.,  directed  energy,  kinetic  energy,  and  neutral  particle  beam 
weapons)  for  defensive  forces.  It  also  has  the  capability  to  model  space-based 
interceptors,  ballistic  missiles,  and  ground-launched  interceptors  as  well  as 
communications,  battle  management,  satellite  attrition  and  sensor  performance.  MEM  is  a 
dynamic,  time-  and  event-step  model.  It  advances  through  events  in  user  specified  time 
increments.  It  is  a  stochastic  model  using  Monte  Carlo  determination  of  results.  MEM 
does  not  require  human  participation  after  the  initial  input  and  output  file  selection. 

[Ref.  7] 

4.  Joint  Integrated  Contingency  Model  (JICM) 

JICM  is  a  theater  level,  multi-sided,  symmetric  and  reactive  model  used  to  support 
balance  assessment,  contingency  analysis  and  is  also  used  as  a  training  tool.  The 
proponent  for  JICM  is  the  Director,  Office  of  the  Secretary  of  Defense  /  Net  Assessment 
(OSD/NA).  It  models  air  and  ground  warfare  in  multiple  theaters,  naval  warfeire,  strategic 
mobility,  joint  and  combined  operations,  and  limited  NBC  warfere.  It  is  a  deterministic 
and  dynamic^  time-step  model.  JICM  can  be  used  as  stand  alone  or  it  can  be  utilized  in  an 
interactive  mode  if  desired.  [Ref.  8] 

5.  Space  Communication  Effectiveness  Model  (SpaceCEM) 

SpaceCEM  is  an  analysis  model  used  to  evaluate  the  performance  of  satellite 

communication  systems.  The  proponent  for  SpaceCEM  is  the  Air  Force  Electronic 
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Warfare  Center,  Special  Studies  Division.  It  models  the  effects  of  jamming,  the 
environment,  terrain  and  nuclear  bursts  on  ground,  air  and  space.  It  models  satellite 
systems  by  a  specific  satellite  or  ground  terminal.  It  is  a  one-sided,  deterministic,  dynamic, 
time-  and  event-step  model.  It  requires  human  decisions  for  processes  and  the  initial  setup 
of  each  anafysis.  However,  it  is  not  an  interactive  model.  [Ref  9] 

6.  Space  Forces  Engagement  Model  (SFEM) 

SFEM  is  an  analysis  model  used  as  a  research  and  evaluation  tool  for  weapon 
systems  development,  systems  effectiveness,  force  capability  and  requirements,  mix 
effectiveness  and  resource  planning.  The  proponent  for  the  SFEM  is  Headquarters,  Air 
Force  Space  Command.  It  models  anti-satellite  weapons  and  maneuver,  decoy,  and 
sabotage  defenses.  It  is  a  two-sided,  asymmetric,  and  dynamic  event-step  model.  It  can 
be  run  either  as  a  stochastic  (Monte  Carlo)  model  or  in  a  deterministic  mode.  Human 
interaction  is  not  allowed.  [Ref  10] 

7.  SIGINT  Analysis  and  Simulation  System  (SASS) 

SASS  is  an  analysis  tool  for  current  and  projected  signal  environments  and  sensor 
systems  operating  in  those  environments.  The  proponent  for  SASS  is  Tactical  Studies  and 
Analysis,  National  Security  Agency.  It  is  a  two-sided,  theater-level  model  that  models 
land,  air,  and  naval  operations.  The  collection  simulation  is  one-sided.  It  can  model 
emissions  Jfrom  the  individual  emitter  level  to  a  mix  of  collection  systems.  It  is  a  dynamic 
time-  and  event-  step  model.  It  models  the  duty  cycle  of  each  emitter  type  stochastically 
using  Monte  Carlo  techniques.  It  requires  human  decisions  for  data  base  ii:q)ut  and 
anatysis  parameters.  [Ref  11] 
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8. 


TACWAR 


TACWAR  is  a  theater  level  model  that  is  used  primarily  for  OPLAN  development 
and  assessment  by  the  CINC  and  by  the  Joint  Staff  for  force  structure  sufficiency 
assessments.  It  is  typically  used  as  a  course  of  action  analyzer  and  a  comparative  analysis 
tool.  It  is  a  deterministic,  two-sided,  dynamic,  time-stepped  model  that  models  air,  land 
and  sea.  The  model  proponent  is  the  Joint  Staff,  Force  Structure,  Resource,  and 
Assessment  Directorate  (J-8),  Automation  Support  Division.  TACWAR  requires  human 
decisions  as  initial  inputs  and  allows  interruption  for  decision  refinements.  [Ref  12] 

9.  THUNDER 

THUNDER  is  a  theater  level  warfare  simulation  model.  It  is  used  as  a  simulation 
to  model  conventional  air,  land  and  naval  air  warfare;  evaluate  force  structures,  conduct 
cost  and  operation  effectiveness  analyses  and  develop  strategies  and  tactics,  among  other 
functions.  The  model  proponent  is  the  Air  Force  Studies  and  Analyses  Agency 
(AFSAA/SAG).  It  is  a  two-sided  model  that  can  be  run  as  a  wargame  or  used  as  an 
analysis  tool.  It  is  a  dynamic,  time-  and  event-  stepped  model,  using  a  Monte  Carlo 
process  to  model  aircraft  kills.  Ground  attrition  is  modeled  deterministically.  No  human 
participation  is  required,  but  are  allowed  if  the  simulation  is  interrupted.  [Ref  13] 

B.  SUMMARY 

From  the  above  and  as  was  concluded  from  Captain  Payne’s  thesis,  there  is  no 
theater  level  simulation  model  that  adequately  models  space.  There  are  space  "tools"  that 
exist,  but  they  are  not  sufficiently  dynamic  to  be  feeder  models  into  theater  and  campaign 
level  models.  This  thesis  will  focus  on  building  a  high  resolution  model  that  can  be 
modified  as  a  feeder  model  for  low  resolution  models  for  analysis  purposes. 
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in.  SIMSATCOM  METHODOLOGY 


The  basic  attributes  of  SIMSATCOM  are  described  in  Section  A.  In  Section  B, 
the  mathematical  formulations,  stochastic  processes,  and  model  products  are  described 
using  the  seven  worksheets  shown  in  Table  1.  A  logic  flowchart  and  detailed  description 
of  the  steps  of  the  computer  simulation,  written  in  Visual  Basic  for  Applications  with 
Microsoft  Excel,  are  provided  in  Appendix  A. 

A.  OVERVIEW 

The  user  initiates  SIMSATCOM  by  defining  the  inputs  for  the  simulation.  These 
input  parameters  are  described  in  Table  1  of  Section  B.  After  the  user  inputs  the 
necessary  parameters,  the  initialcation  portion  of  the  simulation  begins.  Step  1  of  the 
flow  chart  (Appendk  A)  describes  the  procedure  of  how  the  jammer  locations  are 
randomly  generated.  The  user  has  the  option  to  maniially  input  the  jammer  locations  into 
the  “jammer”  sheet  and  turn  off  the  “setjammers”  procedure  in  the  reset  portion  of  the 
program. 

After  the  initialization,  the  new  message  generation  process  begins.  Steps  12 
through  20  in  Appendix  A  describe  how  this  process  operates.  The  simvilation  uses  a 
Poisson  distribution  to  generate  the  number  of  messages  sent  for  each  message  type  for 
each  time  step.  Step  14  describes  how  the  program  uses  the  lambda  values,  input  by  the 
user,  to  generate  the  number  of  messages  for  each  message  type.  Several  Uniform 
random  numbers  are  also  generated  to  create  the  following  data  for  the  new  message: 

•  a  new  location  for  the  sender  and  receiver;  the  sender  and  receiver  type, 
which  indicates  the  type  of  equipment  the  sender  and  receiver  use  to 
counteract  jamming; 
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•  the  number  of  bits,  which  corresponds  to  the  appropriate  type  of  message; 

•  to  determine  if  the  message  is  transmitted  using  the  low  or  medium  data 
rate  if  the  message  is  transmitted  using  satellites  “3”  or  “4”  (models  of 
Milstar  II). 

The  simulation  next  identifies  which  satellite(s)  will  be  transmitting  the  data.  The 
satellite(s)  that  transmit  wUl  depend  on  the  sender  and  receiver  locations.  If  the  sender 
and  receiver  locations  are  not  in  the  same  transmission  area  for  one  satellite  to  transmit, 
then  a  crosslink  is  required,  meaning  that  the  more  than  one  satellite  is  required  to  transmit 
the  message.  Both  satellites  must  be  available  for  the  transmission  to  occur. 

Jamming  of  the  message  may  also  occur  at  both  the  uplink  and  downlink  locations 
of  the  transmission.  There  are  three  possible  outcomes  from  the  jamming  process  in  the 
simulation.  They  are  1)  “totally  jammed”  which  means  that  the  message  was  completely 
jammed;  2)  “partially  jammed”  which  means  that  the  recipient  received  the  message  but 
the  message  was  not  comprehensible;  and  3)  “good”  which  means  that  the  message  was 
not  jammed  and  the  message  was  transmitted  cleanly.  Section  B.5  describes  how  the 
jamming  portion  of  the  simulation  works  in  more  detail. 

If  there  are  available  channels  on  the  appropriate  satellites  and  jamming  does  not 
occur  on  the  transmission,  then  the  message  is  transmitted.  Otherwise  the  message  waits 
in  the  retransmission  queue  for  the  appropriate  time  to  retransmit. 

After  a  message  is  sent,  an  acknowledgment  message  is  transmitted.  This 
acknowledgment  message  informs  the  sender  that  the  message  was  received  successfully 
and  a  retransmission  is  not  necessary.  If  the  acknowledgment  message  is  not  received  by 
the  sender,  then  the  sender  assxunes  that  the  message  was  not  transmitted  successfully  and 
will  retransmit  the  message  again  after  a  set  period  of  time  (an  input  by  the  user,  based 
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upon  the  message  type).  The  sender  will  only  try  to  retransmit  a  message  for  a  set  period 
of  time,  after  which  the  message  transmission  process  stops. 


B.  WORKSHEETS 

Seven  worksheets  are  used  in  the  model.  The  following  section  describes  the 
worksheets  that  are  used  as  data  place  holders.  Table  1  summarizes  each  sheet  function 
and  the  inputs  that  are  necessary  for  initialization. 


Sheet 

Main  Fxmction  Description 

Input  Items 

sim 

Contains  sender  and  receiver  locations, 
satellites  used,  type  of  message  and  status  of 
message 

average  number  of  messages  (lambda  values) 
for  each  message  type 

sat  chan 

Contains  the  maximum  number  of  satellites 
and  the  maximum  number  of  channels  per 
satellite 

number  of  satellites 

maximum  number  of  channels  for  each 
satellite 

time 

Contains  satellite  type  and  corresponding 
baud  rates 

type 

Contains  conversion  chart  for: 

anti-jam  capability 

wait  time  for  sender  to  receive  an 

acknowledgment  message 

maximum  time  the  sender  will  wait  to  send 

a  message 

anti-jam  ability  for  each  ground  station  type 

priority  wait  time  for  each  message  type 

maximum  time  to  wait  to  retransmit  for  each 
message  type 

jammer 

Contains  jammer  locations  and  whether  the 
message  is  jammed 

number  of  jammers 

“perJtot”  value 

jammer  latitude  /  longitude  (optional,  if  input, 
turn  off  reset  jammers  procedure) 

Retrans 

Contains  message  information  similar  to 
‘‘sim”  sheet;  used  as  a  placeholder  for  the 
retransmission  queue 

sat  info 

Contains  satellite  types  and  corresponding 
baud  rates 

satellite  information 

Table  1.  Worksheet  descriptions. 


1.  “Sim” 

The  “sim”  worksheet  (Figure  1)  is  the  principal  worksheet  used.  The  top  portion 
of  the  sheet  contains  input  and  information  placeholder  data.  The  information 
placeholders  include  the  act  rdw  cell  which  shows  that  the  next  available  row  on  the 
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worksheet  is  row  194;  the  cum  time  cell  which  is  the  current  time;  and  the  cum  msgs  cell 
which  denotes  the  number  of  messages  that  have  been  transmitted.  The  input  values  on 
this  sheet  are  the  rate  msg  type  cells  which  are  the  lambda  values.  These  values  are  the 
Poisson  distribution  input  parameters,  which  are  the  average  niunber  of  messages  for  each 
message  type.  These  values  are  used  in  Step  14  of  the  flow  chart  presented  in  Appendix  A 
utilizing  the  Poisson  distribution  shown  in  Equation  1. 


The  next  part  of  the  “sim”  worksheet  contains  the  message  information  used  in  the 
model.  Each  row  represents  a  message  and  the  information  that  accompanies  it.  The  first 
column  represents  the  time  step;  the  second  column  is  the  message  number.  A  negative 
message  number  indicates  that  the  message  is  an  acknowledgment  message  to  the  absolute 
value  of  that  number.  The  message  type  column  indicates  the  type  of  message  that  is 
being  transmitted.  The  possible  values  for  this  column  are  “1”  (strategic  message),  “2” 
(tactical  message),  “3”  (logistics  message),  or  “4”  (administrative  message).  The  method 
used  to  fill  this  cell  is  explained  in  Step  16  of  the  flow  chart  (Appendix  A).  The  number  of 
bits  to  transmit  corre^nds  to  the  type  of  message  transmitting.  The  next  four  columns 
are  the  sender  and  receiver  latitude  and  loi^itude  locations.  The  procedure  for  how  these 
are  iiput  is  in  Step  16  of  the  flow  chart  (Appendk  A). 
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Figure  1.  “sim”  Sheet  Extract 


The  next  two  cells,  the  send  type  and  recv  ’r  type  are  used  for  the  anti-jamming 
capability  of  the  ground  station.  The  value  is  a  random  integer  between  0  and  20.  The 
higher  the  number,  the  better  the  anti-jam  capability  (hence  lower  chance  of  being 
jammed).  The  explanation  of  Figure  4,  the  “type”  worksheet  provides  additional 
information. 

Va&  priority  code  is  currently  set  to  the  same  value  as  the  message  type.  This 
assumes  that  the  smaller  the  message  type  value,  the  higher  the  priority  of  the  message. 
The  random  number  for  the  data  rate  is  used  for  the  data  rate  cell,  which  is  explained  in 
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the  next  paragraph.  The  satellite  type  sender  and  receiver  columns,  sat  type  sender  and 
sat  type  recw,  identify  which  satellites  are  transmitting  the  data.  Table  2  contains  the 
characteristics  and  location  for  each  satellite. 


Satellite  Number 

Satellite  Type 

Longitude 

Possible  Baud  Rate 
{data  rate  column) 

1 

A  (Milstar  I) 

0-90 

LDR 

2 

90  - 180 

LDR 

3 

180-270 

LDR /MDR 

4 

B  (Milstar  II) 

270  -  360 

LDR /MDR 

Table  2.  Satellite  Characteristics. 


To  determine  whether  the  LDR,  “low  data  rate”,  or  MDR,  “medium  data  rate”  is  used,  the 
random  number  generated  earlier  is  compared  to  0.15.  This  value  was  chosen  because  the 
Milstar  II  satellites  have  approximately  15%  MDR  channels  (32  channels)  and  85%  LDR 
channels  (192  channels).  [Ref.  14]  If  the  random  number  is  less  than  0.15  then  MDR  is 
input,  otherwise  the  LDR  is  input. 

The  next  two  columns  specify  data  rates.  The  values  for  these  columns  are  the 
number  of  bits  per  second  that  the  satellite  can  transmit.  The  Up  column  indicates  the 
number  of  bits  transmitting  on  the  uplink.  If  there  is  a  crosslink,  then  the  value  in  the 
crosslink  column  will  be  the  number  of  bits  transmitting  on  both  satellites;  otherwise,  the 
value  in  this  column  will  be  no  crosslink. 

The  TX-baud  value  is  a  function  of  the  bit  rate  and  the  message  size.  It  is  the  total 
number  of  bits  to  transmit  (located  in  the  bits  to  transmit  column  discussed  previously) 
divided  by  the  number  of  bits  per  second  the  satellite  is  able  to  transmit. 

The  jammed  column  indicates  whether  the  message  was  jammed.  The  possible 
values  in  this  column  are  “good”,  indicating  that  the  message  was  transmitted  good,  ‘  j”. 
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indicating  that  the  message  was  not  comprehensible  due  to  jamming,  or  “pj”,  indicating 
that  the  message  was  partially  jammed  (i.e.,  the  message  header  was  comprehensible  but 
the  remaining  portion  of  the  message  was  not).  Further  explanation  is  given  in  the 
“jammer”  sheet  description. 

The  ack  msg  column  is  filled  in  only  if  the  message  is  an  acknowledgment 
message.  The  retry  message  column  indicates  that  the  message  is  a  retransmission.  If  it  is 
a  retransmission,  then  the  cell  is  filled  with  an  “r”.  The  original  message  status  column  is 
only  utilized  in  conjunction  with  an  acknowledgment  message.  The  value  in  this  column 
indicates  how  the  original  message  was  transmitted.  It  contains  either  the  value  “good”  or 
“pj”. 

The  final  colmnn,  msg  sent,  is  the  final  status  of  the  message.  Table  3  contains  the 
possible  values  that  could  be  entered  in  this  column. 


Meaning 

“s” 

message  sent  “good”;  the  sender  received  an  acknowledgment 
message 

“si” 

the  message  was  sent  “good”  but  the  sender  has  not  received  an 
acknowledgment  message 

“s2” 

the  acknowledgment  message  was  sent  (this  value  corresponds  only 
to  an  acknowledgment  message) 

“s3” 

the  message  was  successfiilly  transmitted  but  was  partially  jammed 

‘V’ 

the  message  is  in  a  waiting  status  (Le.,  it  is  currently  being 
processed  or  in  the  retransmission  queue) 

“i” 

the  message  is  jammed 

“stopG” 

The  message  process  is  complete  for  this  message.  The  message 
was  sent  okay,  but  the  sender  never  received  an  acknowledgment 
message  and  the  maximum  time  to  keep  trying  to  send  this  message 
has  been  reached. 

“stopNG” 

The  message  process  is  complete.  In  this  case  the  message  was 
never  received  due  to  jamming,  or  not  enough  channels  were 
available  to  transmit  the  message  and  the  maximum  time  to  keep 
trying  to  send  the  message  was  reached. 

Table  3.  Message  Sent  Column  Possible  Values 
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2.  “sat  chan” 

The  “sat  chan”  worksheet  (Figure  2)  contains  the  number  of  users  that  are 
currently  utilizing  the  specified  satellite.  The  model  maintains  this  sheet  to  ensure  that  the 
satellite  is  not  overloaded.  The  top  portion  of  this  sheet  contains  the  maximum  number  of 
satellites  and  the  maximum  number  of  channels  for  each  satellite.  The  satellite  available 
row  indicates  if  the  satellite  has  any  available  channels  to  transmit  messages.  If  the  value 
is  “False”  then  all  channels  are  being  used  at  that  time.  The  next  row  value  indicates  the 
row  that  is  available  to  input  new  information.  The  current  time  value  indicates  the 
current  model  time. 

The  time  column  shows  the  time  corresponding  to  the  satellite  information  in  the 
adjacent  columns.  The  sail  chan  column  indicates  whether  a  channel  is  being  incremented 
or  decremented  on  the  satellite.  The  cum  sat  column  indicates  the  cumulative  number  of 
channels  being  used  on  the  satellite.  These  two  columns  are  repeated  for  each  satellite  in 
the  model. 
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Figure  2.  ^^sat  chan”  Sheet  Extract 
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3. 


“time” 


The  “time”  worksheet  (Figure  3)  contains  the  message  number  that  currently  is 
being  transmitted  on  a  particular  satellite.  This  sheet  also  keeps  track  of  how  much  time 
that  remains  for  each  transmitting  message.  The  columns  on  the  “time”  sheet  are  broken 
into  categories  by  satellite.  The  first  column  is  the  current  message  transmitting  on  a 
channel  on  the  specified  satellite.  The  time  to  transmit  column  denotes  the  total  amount 
of  time  that  message  will  take  to  transmit.  The  remaining  time  left  on  message  number 
column  indicates  the  time  that  remains  to  transmit  the  message. 
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Figure  3.  “time”  Sheet  Extract 


4.  “type” 

The  “type”  worksheet  (Figure  4)  describes  the  anti-jam  capability  of  the  ground 
station.  The  current  sender  and  receiver  location  and  the  sender  and  receiver  type  are 
fi:om  the  active  row  on  the  “sim”  sheet.  These  values  are  transferred  to  this  worksheet  by 
Step  49  in  the  flow  chart  (Appendix  A). 

Prior  to  the  start  of  the  model,  the  user  inputs  the  proportions  that  correspond  to 
the  type.  The  proportions  represent  the  anti-jamming  capability  of  the  ground  station. 

For  example,  the  sender  in  Figure  4  is  type  8,  which  has  a  0.4  anti-jam  capability.  The 
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probabiKty  that  the  ground  station  will  be  jammed  is  (1  -  0.4)  times  the  probability  of  jam. 
The  ejq)lanation  of  the  “jammer”  sheet  provides  further  details. 

This  sheet  also  contains  the  amount  of  time  the  sender  waits  for  the  receipt  of  an 
acknowledgment  message  prior  to  retransmitting  a  message.  The  values  “1”,  “2”,  “3”, 
and  “4”  correspond  to  the  priority  of  the  message,  and  the  wait  time  for  each  value  is 
input  by  the  user  prior  to  execution  of  the  model.  The  maximum  wait  time  for 
retransmission  is  also  input  for  each  message  priority  class.  This  is  the  maximum  amount 
of  time  the  sender  will  keep  trying  to  send  the  message. 
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5.  “jammer” 

The  “jammer”  worksheet  (Figure  5)  contains  the  number  of  jammers  that  are 
capable  of  jamming  a  particular  site.  This  sheet  is  used  for  every  transmission  to 
determine  if  the  transmission  was  jammed.  The  top  portion  of  the  sheet  contains  the 
current  sender  and  receiver  location  and  the  sender  and  receiver  type  from  the  active  row 
on  the  “sim”  sheet.  These  values  are  updated  by  Step  49  in  the  flow  chart  (Appendbc  A). 
The  value  perJtot  is  used  to  determine  the  portion  of  the  jammed  value  that  results  in  a 
partially  jammed  message,  with  (1- perJtot)  resulting  in  a  totally  jammed  message.  The 
maximum  number  of  jammers  value  is  input  by  the  user  prior  to  the  start  of  the  model 


execution. 


■■■ 

^^■11 

[jSHESISI 

0.7 

Max#  of  jammers  1 

current 

9 

long 

50 

8 

recvr  lat 

38 

Iona 

331 

1 

Jammer  Locations 

U 

plink  1 

Barrier 

lat 

long 

Max  range  of 
jammer 

distance  from 
Jammer  uplink 

probjam 

(uL) 

percjam 

final 

P(jam)  uL 

jam 

none 

jam 

partial 

jam  total 

rand 

draw 

type  jam,  0- 
none, 1-part, 
2-total 

RESULT 

0.99752 

0.6 

0.598513 

0.40149 

0.58104 

1 

0.0711 

0 

1 

-18 

221 

4336 

18340.73617 

0 

2 

-21 

182 

7091 

845.0255738 

0.9858 

3 

-57 

300 

6462 

10058.62184 

0 

4 

-23 

125 

7636 

13570.25114 

0 

5 

-10 

328 

4743 

6 

-5 

8042 

18965.9238 

0 

7 

40 

252 

2808 

1759.273846 

0.60747 

8 

-64 

272 

1948 

5578.250762 

0 

9 

-14 

8 

5816 

3877.799506 

0.55545 

downlink 

overall 

dist  from  jammer 
(downlink) 

prob  jam 
(dL) 

perc 

jam 

final 

PQamJdL 

jam 

none 

jam 

partial 

jam  total 

rand 

draw 

type  jam,  0- 
none,1-part, 
2 -total 

0.95 

0.942752 

0.0572 

0.3401 

1 

2 

0 

10873.55848 

0 

2803.828525 

10667.75358 

954.8089465 

9448.192752 

0 

0 

14810.82914 

14100.92044 

IHHi 

\wmm 

mmm 

F^ure  5.  “jammer”  Sheet  Extract 


19 


In  the  result  portion  of  the  worksheet,  the  top  row  is  a  consolidation  of  the 
jammers  versus  the  sender  and  then  the  receiver.  The  first  column  represents  the  jammer 
nximber.  The  second  and  third  columns  represent  the  location  of  the  jammer.  The  “Reset” 
procedure  runs  the  “setjammers”  routine,  which  assigns  random  locations  for  each  of  the 
jammers.  This  is  explained  further  in  Step  1  of  the  flow  chart  (Appendix  A).  The  next 
column  is  the  maximum  range  of  the  jammer.  All  of  these  values  are  input  at  the  reset 
portion  (Step  1)  of  the  flow  chart.  The  distance  from  the  jammer  uplink  column 
determines  how  far  the  jammer  is  fi:om  the  sender  location.  This  value  is  calculated  by  the 
following  equations:  [Ref  19] 


dist  =  2*  p*  arcsin( 


(2) 


with 

D  =  ^|(x^-X2f  +(j^i -y2f  +(^1  (3) 

Xj  =  /?*sin(^i)*cos(^j) 

=  /?  *  siri(<^, )  *  sin(^i )  (4) 

Zj  =  p*cos(<^,) 


X2=  p*  sin(^2  )  *  cos(02 ) 

y2=  p*sm(^2)*  sin(^2 )  (3) 

Z2  =p*cos(^2) 


j90-latitudel  if  latitude\>0 
*  [90  + latitude!  otherwise 

^0- latitude!  if  latitude2>0 
^  [90+ latitude!  otherwise 
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=36^— longitude!  (8) 

02  =  360  -  longitude!  (9) 

and  p  =  6367  km  (the  radius  of  the  earth). 

If  the  distance  of  the  jammer  to  the  sender  is  less  than  the  maximum  range  of  the 
jammer,  the  probability  of  jam  for  each  individual  location  that  is  within  range  of  the 
jammer  is  determined  by  the  following  formula  which  is  from  the  ACES  model  [Ref  3]: 

Pindiv.jammerCjam)  =  1 -(dist/max  Rauge)^  (10) 

Otherwise  the  probability  of  jam  equals  0.  The  overall  probability  of  jam  is  output  in  the 
Result  row.  The  value  of  the  initial  probability  of  jam  is: 

Pimt(jam)  =  1-(1-Pjam  (first  janimer))*(l-Pjam  (second  jammer))*. ..*(1-Pjam(nth  jammer))  (1 1) 

The  perc  jam  value  is  based  upon  the  type  of  anti-jam  capability  of  the  sender.  The  anti¬ 
jam  effectiveness,  which  is  used  to  compute  the  perc  jam  value,  is  computed  and 
explained  in  Step  50  of  the  flow  chart  (Appendix  A).  The  value  of perc  jam  is: 

perc  jam  =  1  -  anti-jam  effectiveness  (12) 

The  final  probability  of  jam  is  now  calculated  as: 

Pfinai(jam)  =  (PinitOam))  *  perc  jam  (13) 

The  jam  none,  jam  partial,  and  jam  total  values  are: 

jam  none  =  1  -  {jam) 

jam  partial  =  Pfmai  ijarn)  *  {perJtot)  (14) 

jam  total  =  (Jam)  *  (1  -  perJtot) 
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A  random  number  is  then  generated,  as  indicated  in  Step  50  of  the  flow  chart,  to 
determine  whether  the  message  is  totally,  partially  or  not  jammed.  If  the  random  number 
is  less  than  the  jam  none  value  then  the  transmission  is  not  jammed  and  a  “0”  is  entered;  if 
the  random  number  is  between  the  jam  none  and  the  jam  partial  value  then  the 
transmission  is  partially  jammed  and  a  “1”  is  input.  If  the  random  number  is  greater  than 
the  yaw  partial  value,  then  the  transmission  has  been  completely  jammed  and  a  “2”  is 
input. 

For  the  example,  the  sender  is  type  8,  which  corresponds  to  an  anti-jam 
effectiveness  of  0.4.  The  perc  jam  is  0.6  and  the  final  probability  of  jam  is  0.598 
(0.997  *  0.6).  The  random  number  that  was  generated  is  0.071 1.  This  value  is  less  than 
the  jam  none  value  (0.40149),  so  a  “0”  is  entered  in  the  type  jam  ceU. 

Because  both  the  sender  and  receiver  ground  stations  are  susceptible  to  jamming,  a 
similar  calculation  is  needed  for  the  downlink  portion  of  the  transmissioa  The 
confutations  for  the  downlink  columns  work  similarly  to  the  uplink  columns.  However, 
one  difference  is  the  calculation  of  the  type  jam  column.  This  column  is  dependent  upon 
the  result  of  the  uplink  type  jam  column.  If  the  result  from  the  uplink  type  jam  column  is 
jam  total,  then  the  overall  value  for  the  type  jam  column  is  also  yaw  total.  If  the  uplink  is 
partially  jammed,  the  random  number  generated  is  only  conf  ared  with  the  yam  none 
value.  If  the  random  number  is  less  than  the  jam  none  value,  the  value  remains  partially 
jammed  (no  jamming  occurred  on  the  downlink,  but  the  message  was  partially  jammed  on 
the  uplink  so  the  overall  evaluation  is  partially  jammed);  otherwise,  the  message  is  totally 
jammed.  If  the  uplink  was  not  jammed,  then  the  calculation  for  to  colimm  is  the  same 
that  was  used  for  the  uplink  computation  for  fyp§  j^. 

IT 


6. 


“Retrans” 


The  “Retrans”  worksheet  (Figure  6)  contains  the  transmission  information  for  the 
retransmission  queue.  The  messages  that  are  waiting  for  a  retransmission  are  recorded  on 
this  sheet.  The  time  column  indicates  the  time  the  message  was  originally  initiated.  The 
message  number  column  indicates  the  message  number  in  the  queue.  The  message  type, 
bits  to  transmit,  send  lat  and  long,  receiver  lat  and  long,  send  and  receiver  type,  and 
priority  code  columns  are  all  copied  from  the  “sim”  sheet  when  the  message  gets  initiated. 
Step  57  in  the  flow  chart  describes  how  these  values  are  transferred  from  the  “sim”  sheet. 
The  wait  time  value  is  input  based  upon  the  priority  of  the  message,  as  described  in 
Step  56  of  the  flow  chart  (Appendbc  A). 
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Figure  6.  "Retrans”  Sheet  Extract 


The  time  to  try  again  value  is  the  sum  of  the  wait  time  and  the  time  last  try  msg 
column.  The  status  column  indicates  the  status  of  the  retransmission.  “stopG”,  “stopNG” 
and  “good”  values  indicate  that  the  message  will  no  longer  be  retransmitted.  The  value 
“stopG”  means  that  the  original  message  was  sent  successfully  but  the  sender  never 
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received  an  acknowledgment  message.  The  “stopNG”  value  means  that  the  message  was 
never  sent  due  to  jamming  or  no  channel  was  available.  The  value  “good”  indicates  that 
the  message  and  an  acknowledgment  message  were  transmitted  successfully.  The  value 
“w”  indicates  the  message  is  waiting  to  be  retransmitted.  The  value  “n”  indicates  that  the 
message  will  not  be  retransmitted,  (Le.,  the  message  is  a  retry  message  and  the  original 
message  is  already  in  the  queue).  The  time  last  try  message  column  indicates  the  last  time 
the  message  was  atten^ted.  The  satellite  number  up  and  down  columns  indicate  the 
satellite  number  that  will  transmit  the  message.  The  random  number  is  the  data  rate 
random  number  from  the  “sim”  worksheet.  The  jammed  column  indicates  the  original 
message  jam  status.  The  maximum  time  to  try  coliunn  is  the  maximum  amount  of  time  the 
sender  will  try  to  send  the  message.  This  is  updated  as  indicated  by  Step  56  in  the  flow 
chart  (Appendix  A).  The  time  to  stop  trying  is  the  sum  of  the  initial  attempt  to  send  the 
message  {time  column)  plus  the  maximum  time  to  try  the  message.  This  value  indicates 
when  the  model  will  stop  trying  to  send  the  message. 

7.  “sat  info” 

The  “sat  info”  worksheet  (Figure  7)  contains  the  number  and  characteristics  of 
satellites  that  are  used.  The  satellites  are  modeled  after  the  Milstar  I  and  11  satellites.  The 
first  column  indicates  that  the  Type  “A”  satellites  (modeled  after  the  Milstar  I  satellite)  are 
“1”,  “2”,  “5”,  and  “6”  and  Type  “B”  satellites  (modeled  after  the  Milstar  n  satellite)  are 
“3”,  “4”  “7”,  and  “8”.  The  second  column  indicates  the  satellite  number.  The  data  rate 
information  shown  in  the  next  two  columns  gives  the  range  for  the  low  data  rate,  LDR,  as 
75  bits  per  second  to  2,400  bits  per  second.  The  average  is  1237.5  bits  per  second.  For 
MDR,  the  range  for  transmission  is  from  4,800  to  1,544,000  bits  per  second.  The  average 
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rate  is  774,400  bits  per  second.  The  average  rates  are  the  values  used  in  the  calculation  on 
the  “sim”  sheet  for  the  time  to  transmit  a  message. 
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Figure  7.  “sat  info”  Sheet  Extract 


C.  SUMMARY 

After  the  program  has  run,  statistical  procedures  were  used  to  collect  the  data  for 
analysis.  An  example  run  is  shown  in  Chapter  IV  to  show  how  the  model  works,  and  an 
example  of  how  the  model  can  be  used  for  analysis  is  shown  in  Chapter  V. 
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IV.  DEMONSTRATION  CASE 


The  Demonstration  Case  of  the  model  is  an  initial  verification  run  to  ensure  that  all 
of  the  attributes  in  the  model  work  as  they  were  designed.  Four  messages  are  extracted 
firom  the  Demonstration  Case  run  to  exhibit  the  various  features.  This  chapter  describes 
the  features  of  the  model  in  Section  A,  with  the  paths  of  the  messages  detailed  in  Section 
B.  The  paths  of  these  messages  explicitly  describe  how  the  model  works. 


A.  FEATURES 

There  are  several  features  in  the  model  that  are  depicted  as  follows: 

•  Messages  using  an  uplink  only. 

•  Messages  using  a  crosslink,  i.e.,  an  uplink  and  a  downlink  to  transmit  the 
message. 

•  Combinations  of  the  original  message  and  the  acknowledgment  message 
being  transmitted  “good”,  meaning  no  jamming  interfered  with  the 
transmission;  partially  jammed,  meaning  some  of  the  transmission  was 
being  jammed;  or  totally  jammed,  meaning  that  none  of  the  message  was 
received.  For  the  partially  jammed  message,  the  assumption  is  that  only  the 
message  header  was  able  to  transmit  cleanly,  meaning  that  the  receiver 
obtained  the  message,  but  it  was  not  conprehensible. 

•  The  maximum  time  the  sender  waits  to  send  a  message.  This  occurs  if  a 
message  has  been  on  the  retransmission  queue  for  a  specified  amount  of 
time  and  has  not  been  transmitted.  In  this  case  the  sender  stops  trying  to 
resend  the  message  and  it  is  automatically  taken  off  the  queue. 

•  The  maximum  channel  capacity  of  the  model. 
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B.  THE  RUN 


Messages  25, 27,  28,  and  30  are  tracked  for  the  initial  run.  The  features  that  each 
message  cover  in  the  model  are  shown  in  Table  4.  As  the  run  progresses,  the  entries  will 
be  updated. 


Message 

uplink 

only 

uL  & 
dL 

original  msg 
status 

retry  msg 

1  status 

no  avail 
chan,  NAC 

max  time 
expired 

25 

X 

j 

j 

X 

27 

X 

j 

J _ 

X 

28 

X 

NAC 

good 

X 

30 

X 

good 

good 

Table  4.  Summary  of  Features. 


1.  Input 

Inputs  to  the  Demonstration  Case  run  are  shown  in  Table  5. 


srm  sheet 

sat  chan  sheet 

type  sheet 

jammer  sheet 

program  inputs 

lambda  1=  2 

max  num  of 
satellites  =  4 

proportion  of  anti-jam 
capability 

from  type  1  to  type  20  start 
at  0.05  and  increase  by  0.95 
for  each  type 

max  num  of 
jammers  =  9 

num  of 

iterations  =15 

lambda  2  =  3 

maximum  number 
of  channels  on  each 
satellite  =  8 

type  of  msg 

wait  time  /  max  wait  time 

bits  for  type  of 
message: 

lambda  3=1 

1 

6/11 

strategic  = 
2,0OCM,OOO 

lambda  4  =  3 

2 

10/30 

tactical  = 
2,000-8,000 

3 

15/90 

logistic  = 
10,000-30,000 

4 

50/200 

admin  = 
20,000-50,000 
ack  =  2,000 

Table  5.  Inputs  for  Demonstration  Case  Run. 


2.  Time  Step  4 

Figure  8  shows  an  extract  of  the  “sim”  sheet  from  iteration  4.  Note  that  at  time 
step  4,  messages  25, 27, 28  and  30  are  entered  into  the  system.  Message  25  is  jammed, 
message  27  is  partially  jammed,  message  30  is  transmitted  good,  and  there  is  not  an 
available  channel  on  satellite  2  for  message  28  to  be  transmitted,  indicated  by  no  entry  in 


28 


the  status  column.  Messages  25, 27  and  28  are  transmitted  on  different  satellites  for  both 
the  up  and  downlink  and  message  30  is  transmitted  using  the  same  satellite  for  both  the  up 
and  downlinks.  Observe  that  in  the  msg  sent  column  the  status  is  “w”  indicating  that  the 
message  is  waiting  to  be  transmitted. 

Figure  9  shows  that  messages  27  and  30  are  transmitting  on  the  appropriate 
satellites.  Figure  10  shows  an  extract  from  the  “Retrans”  sheet  after  iteration  4.  All  four 
messages  have  been  put  on  the  retransmission  queue,  but  Message  28  but  will  wait  only 
one  time  step,  and  then  try  to  retransmit.  Note  that  the  queue  is  sorted  by  the  type  of 
message  and  then  by  time. 
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Figure  10.  “Retrans”  Sheet  Extract  After  Time  Step  4. 

3.  Time  Step  5 


Figure  1 1  shows  an  extract  from  the  “sim”  sheet  after  iteration  5.  Message  28  has 


been  retransmitted  at  time  step  5  and  the  status  is  now  “good”.  Observe  in  Figtire  12  that 
message  28  has  been  added  to  the  “time”  sheet  satellite  transmission.  Notice  in  Figure  13 
the  time  since  last  try  column  for  the  original  message  28  has  been  updated.  Also  note 


that  the  row  corresponding  to  message  28  at  time  step  5  is  not  in  the  queue  due  to  the  n 


in  the  status  column. 
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Figure  11.  “sim”  Sheet  Extract  After  Time  Step  5. 
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Figure  13.  “Retrans”  Sheet  Extract  After  Time  Step  5. 


4.  Time  Step  7 

Figure  14  shows  that  messages  27  and  30  have  been  sent  at  time  step  7.  The 
status  column  now  indicates  “s3”  for  message  27  and  “si”  for  message  30.  Recall  that 
“s3”  indicates  that  the  original  message  was  a  partially  jammed  message,  and  “si” 
indicates  that  the  original  message  was  sent  “good”.  Notice  that  messages  -30  and  -27  are 
being  transmitted  at  time  step  7.  The  negative  value  corresponds  to  the  message  that  is 
being  acknowledged. 

Figure  15  shows  the  “time”  sheet  after  time  step  7.  Messages  27  and  30  have  been 
taken  off  the  transmission  queue.  The  acknowledgment  messages  -27  and  -30  have  been 
added.  Note  that  there  has  been  no  change  on  the  retransmission  queue  (Figure  1 6) 
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because  the  messages  will  not  be  removed  until  the  acknowledgment  messages  have  been 


received. 
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Figure  14*  “slm”  Sheet  Extract  After  Time  Step  7. 
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Figure  15.  ‘‘time”  Sheet  Extract  After  Time  Step  7. 
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Figure  16.  “Retrans”  Sheet  Extract  After  Time  Step  7. 
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5.  Time  Step  9 

Figure  17  shows  an  extract  from  the  “sim”  sheet  after  iteration  9.  Message  28  has 
been  transmitted  “good”  as  indicated  by  the  “si”  in  message  sent  column  of  the  original 
message  28  row.  The  corresponding  acknowledgment  message  has  been  sent,  which  is 
indicated  by  message  -28.  The  acknowledgment  message  -30  was  sent  as  denoted  by  the 
“s2”  in  the  message  sent  column.  The  message  process  for  message  30  is  now  complete 
as  indicated  by  the  “s”  in  the  message  sent  column.  The  message  is  also  taken  off  the 
retransmission  queue  (Figure  19)  as  indicated  by  the  “good”  in  the  status  column.  Notice 
on  the  “time”  sheet  (Figure  18)  that  message  -28  is  added  and  messages  -27  and  -30  are 
removed.  However,  since  message  27  was  partially  jammed,  it  is  not  taken  off  the 
retransmission  queue  (Figure  1 9)  and  the  status  does  not  change  on  the  “sim  sheet 
message  sent  column. 
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Figure  17.  “sim”  Sheet  Extract  After  Time  Step  9. 
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Figure  18.  “time”  Sheet  Extract  After  Time  Step  9. 
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Figure  19.  “Retrans”  Sheet  Extract  After  Time  Step  9. 


6.  Time  Step  10 

Figure  20  shows  an  extract  from  the  “sim”  sheet  after  iteration  10.  Message  25  is 


being  retransmitted  and  is  now  being  transmitted  “good”.  Note  the  “r”  in  the  retry 
column,  mHirating  that  the  message  is  a  retransmission.  Message  25  has  been  updated  in 
the  retransmission  queue  (Figure  22).  The  time  since  last  try  column  now  shows  10  , 
which  indicates  the  last  time  this  message  was  atten5)ted  was  time  step  10. 
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Figure  20.  “sim”  Sheet  Extract  After  Time  Step  10. 


Figure  21.  “time”  Sheet  Extract  After  Time  Step  10. 


Figure  22.  “Retrans”  Sheet  Extract  After  Time  Step  10. 
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7.  Time  Step  11 

The  “sim”  sheet  (Figure  23)  shows  that  message  27  is  being  retransmitted; 
however,  it  is  being  partially  jammed.  The  acknowledgment  message,  message  -28,  has 
been  sent  as  shown  by  the  “s2”  in  the  message  sent  column.  The  message  process  for 
message  28  is  now  complete  as  indicated  by  the  “s”  in  the  message  sent  column 
corresponding  to  the  message  28  row  and  the  “good”  in  the  status  column  of  the 
“Retrans”  sheet  (Figure  25).  From  the  “time”  sheet  (Figure  24),  message  -28  has  been 
removed  and  message  27  has  been  added.  Message  27  has  also  been  updated  in  the  time 


last  try  colmnn  on  the  “Retrans”  sheet  (Figure  25)  with  the  time,  “1 1”. 
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Figure  23.  “sim”  Sheet  Extract  After  Time  Step  11. 
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Figure  25.  “Retrans”  Sheet  Extract  After  Time  Step  11. 


8.  Time  Step  12. 

At  time  step  12,  message  25  has  been  sent.  The  “sim”  sheet  (Figure  26)  shows  an 
“si”  in  the  msg  sent  column  indicating  that  the  message  was  sent  good.  Message  -25,  the 
corresponding  acknowledgment  mess^e,  is  now  being  transmitted  but  it  is  being  jammed. 
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Figure  26.  “sim”  Sheet  Extract  after  Time  Step  12. 


Figure  27.  ‘‘time’’  Sheet  Extract  After  Time  Step  12. 
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Figure  28.  “Retrans”  Sheet  Extract  After  Time  Step  12. 


9.  Time  Step  15. 

The  final  time  step  for  the  Demonstration  Case  run  is  time  step  15.  The  message 
process  for  messages  25  and  27  are  now  complete;  the  maximum  time  that  the  sender  will 
attempt  to  send  these  messages  has  transpired.  The  “stopG”  in  the  message  sent  column 
of  the  “sim”  sheet  (Figure  29)  and  the  status  column  of  the  “Retrans”  sheet  (Figure  31) 
indicate  that  message  25  was  sent  good,  but  the  sender  never  received  an  acknowledgment 
from  the  receiver.  The  “stopNG”  indicates  that  message  27  was  never  obtained  by 
receiver;  in  this  case  it  was  due  to  the  message  being  jammed.  At  this  time  step,  the 
message  process  for  the  four  messages  being  tracked  is  conqilete,  either  because  the 
messages  were  transmitted  good  or  because  time  expired. 


39 


Figure  29.  “sim”  Sheet  Extract  After  Time  Step  15. 
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Figure  31.  “Retrans”  Sheet  Extract  After  Time  Step  15. 


C.  SUMMARY 

Tables  6, 1,  and  8  show  a  summary  of  the  “sim”,  “time”  and  “Retrans”  sheets, 
respectively,  for  message  27  at  times  4, 7,  9,  10, 11, 14,  and  15. 


time 

status 

4 

Message  generation,  jammed  value  =  “pj” 

7 

Afeg  sent  value  changed  to  “s3”,  message  -27  (acknowledgment  message  for  message  27  is  sent) 

IB 

Message  -27  has  transmitted  (“s2”  value  in  msg  sent  col);  Message  27  status  has  not  changed 
since  message  was  partially  jammed  and  it  needs  to  be  retransmitted 

There  were  no  available  channels  (satellites  2  and  3  are  at  capacity)  so  no  retransmission 

11 

Message  retransmitted,  jammed  value  =  “pj”,  retry  from  retrans  column  value  =  “r” 

14 

Message  -27  is  jammed 

15 

Message  process  complete;  msg  sent  value  =  “stopNG”  indicated  message  was  never 
transmitted  successfidly 

Table  6.  Sheet  Summary  for  Message  27. 


time 

status 

4 

Added  to  “time”  sheet  queue,  transmitting  on  satellites  2  and  3 

7 

Message  27  removed  from  “time”  sheet  queue,  message  -27  added 

9 

Message  -27  removed  from  “time”  sheet  queue 

11 

Added  to  “time”  sheet  queue,  transmitting  on  satellites  2  and  3 

14 

Message  27  removed  from  “time”  sheet  queue 

Table  7.  “time”  Sheet  Summary  for  Message  27. 


time 

status 

4 

Added  to  retransmission  queue;  time  to  try  a^ain  value  =  10,  time  to  stop  trying  =15 

11 

Time  to  try  again  value  =  17,  time  last  try  msg  =11 

15 

Time  to  stop  trying  value  =15,  status  =  “stopNG”  a  flag  that  the  message  process  is  complete 
for  this  message. 

Table  8.  “Retrans”  Sheet  Summary  for  Message  27. 
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This  chapter  described  an  example  run  with  example  input  parameters.  There  are 
many  parameters  that  can  be  changed  to  answer  a  variety  of  questions.  Chapter  V  will 
choose  two  such  parameters  to  show  how  the  model  can  be  used  for  data  analysis  to 
answer  specific  questions. 
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V. 


ANALYSIS 


The  purpose  of  this  chapter  is  to  state  the  measures  of  effectiveness,  describe  the 
cases  that  were  analyzed  and  describe  how  the  analysis  was  conducted.  There  are  many 
model  parameters  that  can  be  modified  for  analysis;  however,  the  focus  of  this  thesis  is  the 
examination  of  two  parameters:  anti-jam  capability  and  the  number  of  channels  available 
per  satellite,  reflected  in  the  MOEs  defined  in  Section  A  below. 

Four  cases  are  examined  to  complete  this  analysis  as  follows: 

•  Case  1  (baseline  model  for  the  analysis)  -  the  anti-jam  capability  and 
channel  capacity  are  set  at  low  levels  (LL); 

•  Case  2  -  the  number  of  available  channels  is  doubled  fi'om  the  baseline 
model,  but  the  anti-jam  capability  remains  low  (HL); 

•  Case  3  -  the  anti-jam  capability  is  increased  from  the  baseline  model,  but 
the  channel  capacity  remains  low  (LH); 

•  Case  4  -  the  number  of  channels  is  doubled  and  the  anti-jam  capability  is 
increased  firom  the  baseline  model  (HH). 

Model  run  results  are  described  in  detail  in  Section  B.  A  2  x  2  Analysis  of  Variance 
(ANOVA)  is  used  in  the  analysis  to  compare  these  four  cases.  A  detailed  discussion  of 
the  analysis  methodology  and  results  are  provided  at  Section  C. 

A.  MEASURES  OF  EFFECTIVENESS  (MOE) 

To  help  answer  the  question,  “Would  it  be  better  to  buy  a  satellite  system  that  has 
a  better  anti-jam  capability  or  has  more  channel  capacity?”,  the  following  measures  of 
effectiveness  are  identified: 

•  The  total  number  of  good  messages  transmitted. 
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•  The  number  of  messages  that  were  transmitted  successfully,  but  were 
stopped  because  a  successful  acknowledgment  message  was  not  received 
before  the  maximum  time  to  attempt  transmission  had  e5q)ired  (stopG). 

•  The  number  of  messages  that  were  never  transmitted  (either  due  to 
jamming  or  no  available  channels  when  the  message  process  stopped) 
because  the  maximum  time  to  attempt  transmission  had  expired  (stopNG). 

•  The  average  time  to  send  a  message,  given  that  it  was  transmitted 
successfully. 

•  The  proportion  of  good  messages  to  the  number  of  messages  generated  for 
each  priority. 

B.  CASES 

The  cases  were  selected  to  determine  if  changing  the  model’s  parameters  would 
cause  a  significant  difference  in  the  MOEs  (i.e.,  if  the  anti-jam  capability  was  increased, 
the  number  of  channels  was  increased,  or  both  were  increased).  The  random  number  seed 
was  set  to  a  different  value  at  the  beginning  of  each  run  for  each  case.  Each  case  was  run 
5  times  for  60  simulated  seconds  for  each  run.  For  an  Intel  Pentium  133  MHz  processor, 
each  run  takes  approximately  2  hours  and  30  minutes. 

1.  Case  1  (LL) 

The  input  parameters  for  Case  1  are  given  in  Table  9.  Both  fectors  that  are  being 
conpared  in  Case  1  are  set  at  (low,  low)  (i.e.,  the  channel  capacity  is  set  at  eight  channels 
and  the  anti-jam  capability  is  set  at  0.25). 

Figures  32  and  33  show  the  number  of  messages  generated  by  priority  over  time 
and  the  average  number  of  messages  generated  over  time,  respectively.  These  charts 
show  that  the  number  of  messages  generated  from  the  Poisson  Distribution  for  each  type 
of  message  approach  the  respective  lambda  values  given  in  Table  9.  Appendk  D  (Figures 
97  -  102)  shows  similar  figures  for  Cases  2, 3  and  4  with  similar  results. 
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The  total  number  of  Case  1  messages  that  were  transmitted  successfully,  “stopG” 
and  “stopNG”  over  time  are  shown  in  Figure  34.  The  number  of  channels  used  for  each 
satellite  over  time  is  shown  at  Figure  35,  which  shows  that  almost  aU  of  the  channels  on 
each  satellite  are  being  used.  Similar  graphs  that  show  channel  use  by  priority  are  included 
in  Appendfac  C  (Figures  77-81). 

Figure  36  shows  the  average  time  to  transmit  a  message  for  every  good  message 
transmitted  by  priority.  As  time  increases,  the  time  to  transmit  a  message  depends  upon 
the  priority  of  the  message.  Because  the  size  of  the  message  is  based  on  the  type  of  the 
message,  the  results  correlate  to  the  input  parameters. 

Figure  37  shows  the  proportion  of  good  messages  of  each  priority  to  the  number 
of  messages  generated  for  each  respective  priority  (i.  e.,  number  of  priority  1  good 
messages  /  number  of  total  priority  1  messages  generated).  The  message  with  the  higher 
priority  will  have  a  higher  probability  of  being  transmitted  successfully.  The  output 
correlates  to  the  model  design  which  processes  messages  first  by  priority,  followed  by 
first-in,  first-out. 


“sim”  sheet 

“sat  chan”  sheet 

“type”  sheet 

“jammer” 

sheet 

Other 

— 

lambda  1=2 

lambda  2  =  3 

lambda  3  =  1 
lambda  4  =  3 

Maximum  number 
of  satellites  =  4 
Maximum  number 
of  channels  on  each 
satellite  =  8 

Anti-jam  ability 

Type  of  message 

1 

2 

3 

4 

For  each  type  =  25 

Wait  time  /  max  wait  time 

3/7 

6/13 

9/19 

12/25 

Maxnumof 
jammers  =  9 

, 

Num  of  iterations  =  60 

Bits  for  type  of  message 

Strategic:  1,000  -  2,000 
Tactical:  1,000  -  4,000 
Logistic:  4,000  -  6,000 
Administrative:  4,000  - 
8,000 

Acknowledgment:  2,000 
Random  number  seeds 
Runs  1-5: 

1973272912 

281629770 

20006270 

1280689831 

2096730329 

Table  9.  Baseline  Input  Parameters. 
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No.  of  Messages  |  No.  of 


F^ure  32.  Case  1  (LL):  Number  of  Messages  Generated  by  Priority. 
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Figure  33.  Case  1  (LL):  Average  Number  of  Messages  Generated. 
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@#good  ■#stopG  □#stopNG 


Figure  34.  Case  1  (LL):  Cumulative  Message  Status  Summary, 
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Figure  35.  Case  1  (LL):  Cumulative  Number  of  Channels  in  Use  Over  Time. 
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proportion 


Figure  37.  Cumulative  Proportion  of  Good  Messages  to  Number  of  Messages  Generated  by  Priority. 
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2.  Case  2  (HL) 

The  input  parameters  for  Case  2  that  were  modified  from  Case  1  are  channel 
capacity  and  the  random  number  seeds  for  each  run.  The  factors  are  now  at  (high,  low) 
levels.  The  number  of  available  channels  on  each  satellite  is  increased  from  8  to  16.  The 
random  number  seeds  for  runs  1-5  were:  1933576050, 913566091,  246780520, 
1363774876,  and  604901985.  All  other  input  parameters  remained  the  same  as  those  in 
Case  1. 

The  total  number  of  Case  2  messages  that  were  transmitted  successfully,  “stopG” 
and  “stopNG”  over  time  are  shown  in  Figure  38.  Table  10  shows  a  comparison  of  good 
messages  transmitted  at  different  time  intervals  (the  values  are  extracted  from  the  number 
of  good  messages  shown  in  Figures  34  and  38.)  The  number  of  good  transmissions  from 
Case  2  is  at  least  1 .3  times  greater  than  Case  1  at  the  specified  time  intervals. 

The  number  of  channels  used  for  each  satellite  over  time  is  shown  at  Figure  39. 
Similar  graphs  that  show  channel  use  by  priority  are  included  in  Appendfac  C  (Figures  76  - 
80).  Not  all  channels  of  all  satellites  are  being  used.  This  could  be  due  to  several  fectors: 
1)  at  the  time  messages  were  being  transmitted,  they  were  being  jammed;  2)  one  of  the 
satellites  required  to  be  used  to  transmit  a  message  was  at  full  capacity,  and  the  other 
satellite  was  available;  or  3)  satellites  were  available,  but  messages  were  in  the  wait 
portion  of  the  retransmission  queue.  The  lack  of  channels  at  full  capacity  is  more  apparent 
in  Case  2  than  Case  1  because  each  of  the  fectors  mentioned  above  is  more  pronoxmced  in 
Case  2.  In  particular,  there  are  more  channels  available  in  Case  2  and  therefore,  more 
chances  for  each  of  the  feictors  to  occur. 
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Figure  40  shows  the  average  time  to  transmit  a  message  for  every  good  message 
transmitted  by  priority.  Similar  to  Case  1,  as  time  increases,  the  time  to  transmit  depends 
upon  the  priority  of  the  message.  Because  the  size  of  the  message  is  based  on  the  type  of 
the  message,  the  results  correlate  to  the  input  parameters.  Compared  to  Case  1,  the  time 
to  transmit  priority  1, 2  and  3  messages  went  down  slightly,  but  the  time  to  transmit 
priority  4  messages  increased  slightly  (0.2  seconds). 

The  proportion  of  good  messages  to  the  total  number  of  messages  is  shown  in 
Figure  41.  The  proportion  has  increased  from  Case  1  for  all  priorities. 


Time  (seconds) 

Niunber  of  Good  Message  Transmissions 

Proportion  Increase  in  Number 
of  Good  Transmissions  in  Case  2 
Over  Case  1 

Case  1 

Case  2 

15 

12.8 

17.4 

1.36 

30 

32.8 

54.2 

1.65 

45 

54.6 

89.8 

1.64 

60 

76.2 

125.2 

1.64 

Table  10.  Good  Message  Transmission  Time  Comparison  for  Cases  1  (LL)  and  2  (HL). 
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Figure  39.  Case  2  (HL):  Cumulative  Number  of  Channels  in  Use  Over  Time. 


Figure  40.  Case  2  (HL);  Average  Time  to  Transmit  Good  Messages. 
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Figure  41.  Cumulative  Proportion  of  Good  Messages  to  Number  of  Messages  Generated  by  Priority. 


3.  Case  3  (LH) 

The  input  parameters  that  were  changed  from  Case  1  for  Case  3  are  the  anti-jam 
capability  and  the  random  number  seeds  for  each  run.  The  factors  are  now  at  (low,  high) 
levels.  The  satellite  capacity  is  eight  channels  (low),  the  anti-jam  ability  for  each  type  was 
changed  from  0.25  to  0.75  (high)  and  the  random  number  seeds  for  runs  1-5  were: 
1511192140, 1259851944,  824064364, 150493284,  and  24270853 1 .  All  other  input 
parameters  remained  the  same  as  Case  1. 

The  total  number  of  Case  3  messages  that  were  transmitted  successfully,  “stopG” 
and  “stopNG”  over  time  are  shown  in  Figure  42.  Table  1 1  shows  the  values  for  the 
number  of  good  messages  transmitted  at  15  second  increments.  The  number  of  good 
messages  transmitted  has  increased  over  2.5  times  from  Case  1  (after  6  initialization 
seconds). 
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The  number  of  channels  used  for  each  satellite  over  time  is  shown  at  Figure  43. 
Similar  graphs  that  show  channel  use  by  priority  are  included  in  Appendix  C  (Figures  87  - 
91).  Figure  43  shows  an  approximately  horizontal  line  for  aU  satellites,  which  indicates 
that  almost  all  of  the  channels  are  in  use  as  in  Case  1 . 

Figure  44  shows  the  average  time  to  transmit  a  message  for  every  good  message 
transmitted  by  priority.  Similar  to  Case  1,  as  time  increases,  the  time  to  transmit  depends 
upon  the  priority  of  the  message.  Because  the  size  of  the  message  is  based  on  the  type  of 
the  message,  the  results  correlate  to  the  input  parameters.  The  time  to  transmit  priority  1, 
2  and  4  messages  increased,  and  the  time  to  transmit  priority  3  messages  decreased.  The 
increased  time  to  transmit  messages  may  be  due  to  more  messages  in  the  retransmission 
queue  that  will  be  sent  successfully,  when  the  opportunity  for  transmission  occurs.  The 
analysis  of  variance  discussed  later  shows  that  there  is  no  statistical  significance  to  these 
results. 

The  proportion  of  good  messages  to  the  number  of  messages  generated  for  each 
priority  is  shown  in  Figtire  45.  The  proportion  has  increased  from  Case  1  for  all  priorities. 
The  proportion  of  priority  1  and  2  good  messages  transmitted  rose  substantialfy  compared 
to  Case  1.  The  proportion  of  priority  3  and  4  good  messages  transmitted  did  not  mcrease 
substantially  con:5)ared  to  Case  1.  This  is  because  the  maximum  time  expired  before  there 
was  a  channel  available  to  transmit  them  (see  Appendix  B;  Figures  69  -  70  did  not  change 
substantially  conpared  to  Figures  59  -  60). 
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time 


I  #  good  ■  #  stopG  □  #stopNG 


Figure  42.  Case  3  (LH):  Cumulative  Message  Status  Summary. 


Time  (seconds) 

Number  of  Good  Message  Transmissions 

Proportion  Increase  in  Number  of 
Good  Transmissions  in  Case  3  Over 
Case  1 

Case  1 

Case  3 

15 

12.8 

33.2 

2.6 

30 

32.8 

84.2 

2.8 

45 

54.6 

136.2 

2.5 

60 

76.2 

2.5 

Table  11.  Good  Message  Transmission  Time  Comparison  for  Cases  1  (LL)  and  3  (LH). 
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Figure  45.  Cumulative  Proportion  of  Good  Messages  to  Number  of  Messages  Generated  by  Priority. 


4.  Case4(HH) 

The  input  parameters  that  were  changed  for  Case  4  (HH)  from  Case  1  are 
channel  capacity,  anti-jam  capability  and  the  random  number  seeds  for  each  run.  The 
factors  are  now  at  (high,  high)  levels.  The  number  of  available  channels  on  each  satellite  is 
increased  from  8  to  16  (high);  the  anti-jam  capability  for  each  type  was  changed  from  0.25 
to  0.75  (high);  and  the  random  number  seeds  for  runs  1-5  were:  75253171, 1964472944, 
1202299975, 233217322,  and  191 1216000.  All  other  input  parameters  remained  the 
same  as  Case  1 . 

The  average  number  of  messages  for  all  runs  of  Case  4  that  were  transmitted 
successfully,  “stopG”  and  “stopNG”  over  time  are  shown  in  Figure  46.  Table  12  shows 
the  values  for  the  number  of  good  messages  transmitted  at  15  second  increments.  The 
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number  of  good  messages  transmitted  has  increased  over  4  times  from  Case  1 .  The 
number  of  channels  used  for  each  satellite  over  time  is  shown  at  Figure  47.  Similar  graphs 
that  show  channel  use  by  priority  are  included  in  Appendix  C  (Figures  92  -  96). 

Figure  48  shows  the  average  time  to  transmit  a  message  for  every  good  message 
transmitted  by  priority.  Similar  to  Case  1,  as  time  increases,  the  time  to  transmit  depends 
upon  the  priority  of  the  message.  Because  the  size  of  the  message  is  based  on  the  type  of 
the  message,  the  results  correlate  to  the  input  parameters.  The  time  to  transmit  priority  1 
and  4  messages  increased  from  Case  1,  and  the  time  to  transmit  priority  2  and  3  messages 
decreased. 

The  proportion  of  good  messages  to  the  number  of  messages  generated  for  each 
priority  is  shown  in  Figiire  49.  The  proportion  has  increased  from  Case  1  for  all  priorities. 
Priority  4  messages  did  not  increase  as  substantially  as  priority  1, 2  and  3  good  messages. 
This  is  because  maximum  time  expired  before  there  was  a  channel  available  to  transmit 
them  (see  Appendbc  B;  Figure  75  did  not  change  substantially  conqjared  to  Figure  60). 
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I  #  good  ■  #  stopG  □  #stopNG 


Figure  46.  Case  4  (HH):  Message  Status  Summary. 


Time  (seconds) 

Number  of  Good  Message  Transmissions 

Proportion  Increase  in  Number 
of  Good  Transmissions  in  Case  4 
Over  Case  1 

Case  1 

Case  4 

15 

12.8 

56.4 

4.4 

30 

32.8 

142.2 

4.3 

45 

54.6 

225.8 

4.1 

60 

76.2 

314.8 

4.1 

Table  12.  Good  Message  Transmission  Time  Comparison  for  Cases  1  (LL)  and  4  (HH). 
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No.  of  channels  in  use 


— ♦—Priority  1  Priority  2  Priority  3  Priority  4 

Figure  49.  Proportion  of  Good  Messages  to  Number  of  Messages  Generated  by  Priority. 


5.  Summary 

Stimmary  statistics  for  each  case  and  MOE  were  collected  and  are  shown  at 
Appendix  E.  Table  13  shows  a  summary  for  the  total  number  of  good  transmissions  for 
each  case  over  15  simulated  second  intervals.  The  number  of  good  transmissions  increase 


from  Case  1  based  upon  the  change  of  parameters. 


Time  (seconds) 

Number  of  Good  Messages  Transmitted 

Case  1  (LL) 

Case  2  (HL) 

Case  3  (LH) 

Case  4  (HH) 

15 

12.8 

17.4 

33.2 

56.4 

30 

32.8 

54.2 

84.2 

142.2 

45 

54.6 

89.8 

136.2 

225.8 

60 

76.2 

125.2 

190.4 

314.8 

Table  13.  Total  Number  of  Good  Message  Transmissions  Over  Time. 


Table  14  shows  a  summary  for  the  average  time  to  transmit  messages  for  each  case 
by  priority  at  time  60.  The  results  do  not  show  a  significant  difference  between  the 
different  cases.  The  analysis  of  variance  discussion  in  Section  C  considers  whether  there  is 


a  significant  difference  for  the  average  time  to  transmit. 
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Case  1  (LL) 

Case  2  (HL) 

Case  3  (LH) 

Priorityl 

5.71 

5.47 

6.36 

5.89 

Priority  2 

8.36 

7.96 

7.41 

Priority  3 

12.65 

10.26 

11.98 

Priority  4 

13.45 

13.62 

14.14 

14.43 

All  messages 

9.0 

9.22 

9.09 

8.76 

Table  14.  Average  time  to  transmit  Messages  by  Priority,  Time  60. 


Table  15  shows  a  summary  of  the  proportion  of  good  messages  to  the  number  of 


messages  generated  for  each  case  by  priority  at  time  60.  The  proportion  comparisons 
from  Case  1  to  the  other  three  cases  are  similar  to  the  results  of  those  in  Table  13. 


Case  1  (LL) 

Priorityl 

0.19 

0.26 

0.72 

0.18 

0.24 

0.42 

EIQHHHHH 

ISSBISSHliHii 

0.17 

0.22 

0.14 

OSSHHHHli 

Priority  4 

ebbhhhh 

0.14 

0.35 

0.57 

Table  15.  Proportion  of  Good  Messages  to  Number  of  Messages  Transmitted,  Time  60. 


C.  ANALYSIS  OF  VARIANCE  (ANOVA) 

A  two-way  Analysis  of  Variance  (ANOVA)  was  conducted  for  each  MOE.  The 
following  model  was  used  to  represent  this  two-fector  experiment; 

Xgk  =  +  Pj+Yij+Syk 

where  p  represents  the  expected  value  over  all  popiilations,  a  represents  the  anti-jam 
fector,  P  represents  channel  capacity  fector,  y  represents  the  interaction  parameter 
between  the  two  fectors,  and  s  represents  the  random  error  term.  Factor  a  has  two  levels, 
i  =  1,  2  (for  low  and  high  anti-jam  ability).  Factor  p  has  two  levels,  j  =  1, 2  (for  low  and 
high  number  of  channels);  and  k  =  1, 2, 3, 4,  and  5  observations  [Ref.  16,  p.  415]. 

Table  16  shows  the  ANOVA  for  the  number  of  good  transmissions  for  all 
messages  (a  =  0.05).  The ‘T”  value  is  greater  than  “F  statistic”  value  for  both  fectors 
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indicating  that  there  is  a  difference  between  the  means  for  anti-jam  capability  and  channel 
capacity.  These  results  imply  that  the  main  effects  of  anti-jam  capability  and  the  channel 
capacity  do  affect  the  number  of  good  messages  transmitted.  The  “F”  value  for  the 
interaction  effect  is  greater  than  the  “F  statistic”,  which  implies  that  there  is  also  an 
interaction  effect  between  the  two  fectors.  The  presence  of  an  interaction  effect  is  also 
shown  in  Figure  56.  Appendix  F  shows  the  remaining  ANOVA  results  and  the  interaction 
plots  for  the  MOEs  for  all  messages  (Table  18,  Figures  103,  108, 1 13, 118,  and  123)  and 
each  priority  (Tables  19  -  22,  Figures  104  -  107, 109  -  1 12, 1 14  -  117, 1 19  -  122, 124  - 
127). 


ANOVA  results  for  the  time  to  transmit  for  all  messages  indicate  the  anti-jam 
capability  and  channel  capacity  factors  do  not  affect  the  time  to  transmit  messages 
(Appendix  F,  Table  18  and  Figure  123).  The  ANOVA  results  for  number  of  “stopG” 
messages  (Appendix  F,  Table  1 8)  show  that  there  is  a  difference  between  the  means  for 
both  factors,  but  there  is  no  interaction  effect.  The  parallel  lines  in  the  “stopG”  interaction 


plot  (Appendix  F,  Figure  108)  also  shows  the  no  interaction  effect. 


ANOVA 

Source  of  Variation 

Df 

SS 

MS 

F 

Fcrit 

J  (Anti- Jam) 

1 

115368.1 

115368.05 

265.9169 

4.494 

C  (Channel 

1 

37584.45 

37584.45 

86.63006 

4.494 

Capacity) 

JxC  (Interaction) 

1 

7106.45 

7106.45 

16.37997 

4.494 

Error 

16 

6941.6 

433.85 

Total 

19 

167000.6 

Table  16.  ANOVA  Good  Transmissions  (All  messages). 
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The  ANOVA  results  (Appendix  F,  Figure  1 13)  for  the  “stopNG”  messages 
indicate  that  there  is  a  difference  between  the  means  for  anti-jam  capability  and  channel 
capacity,  and  there  is  an  interaction  effect.  The  ANOVA  results  (Appendix  F,  Figure  118) 
for  the  proportion  of  good  messages  to  the  number  of  messages  generated  for  each 
priority  indicate  that  there  is  a  difference  between  the  means  and  there  is  an  interaction 
effect. 

Table  17  displays  a  summary  of  the  ANOVA  results  for  each  MOE  for  all 
messages. 

The  Analysis  of  Variance  is  a  preliminary  data  analysis  result  for  this  model.  The 
focus  of  this  thesis  is  not  to  conduct  data  analysis  using  the  model,  but  to  describe  and 
demonstrate  the  features  of  the  model  and  how  it  can  be  used  to  answer  specific  questions. 
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MOE 

Difference  Between  Means 

Interaction  Effect 

Number  of  good 
transmissions 

YES 

YES 

Number  of  “StopG” 
messages 

YES 

NO 

Number  of  “StopNG” 
messages 

YES 

YES 

Time  to  transmit  messages 

NO 

NO 

Number  of  good  messages  / 
number  of  total  messages 

YES 

YES 

Table  17.  ANOVA  Results  Summary  (For  all  messages). 
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VI.  RECOMMENDATIONS 


The  methodology  and  model  design  presented  in  this  thesis  provide  a  useful  model 
for  answering  questions  about  communications  involving  geosynchronous  satellites.  Many 
features  of  the  model  are  useful  for  answering  a  variety  of  questions.  However,  there  are 
some  areas  of  the  model  where  it  would  be  beneficial  to  add  additional  features  to  provide 
further  analysis  capabilities. 

The  model  is  not  dynamic  in  the  area  of  satellite  locations.  It  would  be  beneficial 
to  modify  the  model  in  order  for  the  user  to  input  the  location  of  the  satellites  and  to 
designate  specific  areas  of  responsibility  for  each  satellite.  In  building  this  modification,  it 
would  be  useful  to  have  satellites  with  different  capabilities  at  specified  locations  (i.e., 
different  data  rates,  ranging  jfrom  low  to  high,  modify  the  anti-jam  capability  for  each 
satellite  and  have  specified  fi’equencies  on  each  satellite). 

It  would  also  be  advantageous  to  modify  a  portion  of  the  input  locations  to  a 
specific  area  of  interest.  This  could  be  done  with  a  slight  modification  of  the  current 
model  (i.e.,  add  a  percentage  that  would  set  to  a  specified  area  of  interest  and  the 
remainder  to  the  upper  and  lower  bounds  for  the  input  latitude  and  longitude). 

A  more  sophisticated  jamming  process  could  also  be  added  to  the  model. 

Different  types  of  jammers  with  different  capabilities  would  be  a  beneficial  improvement. 
For  example,  a  microwave  radar  can  cause  interference  to  communication.  A  simulated 
microwave  radar  could  be  a  type  of  jammer  with  appropriate  firequencies  that  it  could  jam. 
After  adding  more  sophisticated  jammers,  it  would  be  usefial  to  allow  the  sender  and 
receiver  to  select  and  change  fi-equencies  if  a  specific  frequency  is  jammed. 
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Once  these  additions  are  made,  the  simulation  could  be  added  as  a  high  resolution 
module  to  other  models.  To  make  this  model  computable  with  another  model,  the  input 
portion  of  the  model  would  require  modification.  The  jammers  and  the  output  portion  of 
the  model  woxdd  also  need  to  be  tailored  to  the  parent  model. 

Models  that  require  near-earth  or  semi-synchronous  satellites  could  be  designed  in 
a  similar  fashion.  However,  the  location  of  the  satellites  would  be  dependent  on  the 
location  of  the  satellite  at  the  start  of  the  simulation  and  the  time  in  the  model  a  particular 
action  needs  to  take  place.  Specifications  specific  to  the  satellites  such  as  sensor 
capabilities  and  their  corresponding  equations  would  need  to  be  added. 

Finally,  JWARS  would  be  able  to  utilize  this  model  to  investigate  the  satellite 
aggregation  issue.  With  some  modification,  it  could  become  the  high  resolution  module 
for  satellite  representation  in  JWARS. 
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APPENDIX  A.  SIMSATCOM  PROGRAM  DESCRIPTION 

This  appendix  discusses  random  number  generation,  flow  chart  model  logic  and  a 

detailed  description  of  each  step.  The  flowcharts  shown  in  Figures  51-56  describe  the 

steps  of  the  simulation.  An  extract  of  the  Uniform  Random  Number  Generator  method 

follows.  The  method  is  converted  from  the  Pascal  Random  Number  Generator  described 

in  Law  and  Kelton.  [Ref  18,  p.  451-454] 

Zi  =  stream 
Hil5  =  Zi\B2E16 

Lowprd  =  (Zi  -  Hil5  *  B2E16)  *  MULTI 
Lowl  5  =  Lowprd  \  B2E 1 6 
Hi31  =  Hil5  ♦  MULTI  +  LowlS 
Ovflow  =  Hi31\B2E15 

Zi  =  (((Lowprd  -  LowlS  *  B2E16)  -  MODLUS)  +  (Hi3 1  -  Ovflow  *  B2E15)  *  B2E16)  +  Ovflow 
If  Zi  <  0  Then  Zi  =  Zi  +  MODLUS 
Hil5  =  Zi\B2E16 

Lowprd  =  (Zi  -  HilS  *  B2E16)  *  MULT2 
LowlS  =  Lowprd  \  B2E16 
Hi31  =  Mils  *  MULT2  +  LowlS 
Ovflow  =  Hi31\B2ElS 

Zi  =  (((Lowprd  -  LowlS  *  B2E16)  -  MODLUS)  +  (Hi31  -  Ovflow  ♦  B2E1S)  ♦  B2E16)  +  Ovflow 
If  Zi  <  0  Then  Zi  =  Zi  +  MODLUS 
stream  =  Zi 

mdm  =  (2  *  (Zi  \  2S6)  +  1)  / 16777216 
The  following  subsections  describe  the  steps  of  the  simulation. 

1.  Reset 

The  “reset”  procedure  sets  the  sheets  on  the  spreadsheet  to  the  starting  position. 
It  also  runs  the  procedures  to  set  the  number  of  satellites  and  the  number  of  channels  on 
each  of  the  satellites.  This  updates  the  entire  “sat  chan”  sheet  with  the  information  that 
was  input  by  the  user.  “Reset”  also  runs  the  “setjammers”  routine,  which  assigns  random 
locations  for  each  of  the  jammers.  An  extract  of  the  “setjammers”  procedure  follows: 


For  numJam  =  1  To  maxNumJam 
UhifRand  mdNum,  stream 

Cells(miim,  latcol).  Value  =  Int((upbndLat  -  lowbndLat  +  1)  *  mdNum  +  lowbndLat) 
UnifEland  mdNum,  stream 
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Cells(mum,  longCol). Value  =  Int((upbndLong  -  lowbndLong  +  1)  *  mdNum  + 
lowbndLong) 

UnifRand  mdNum,  stream 

Cells(mum,  maxmgcol). Value  =  Int((upbndRng  -  lowbndRng  +  1)  *  mdNum  + 
lowbndRng) 
mum  =  mum  +  1 

Next 

If  the  user  desires  to  employ  fixed  locations  for  the  jammers,  this  procedure  in  the 
“reset”  routine  must  be  commented  out  and  the  locations  manually  entered  on  the 
“jammer”  sheet  prior  to  running  the  program. 

2.  Execute  the  New  Message  Process 

Step  2  initiates  the  new  message  process  at  Step  12.  After  the  execution  of  this 
procedure,  continue  with  Step  3. 

3.  Number  of  Iterations 

This  is  a  counter  for  a  ‘Tor”  loop  to  run  the  model  for  the  set  number  of  time 
increments,  with  each  increment  being  one  second. 

4.  Increment  Time 

Once  all  new  messages  have  been  processed,  time  on  the  “sim”  sheet  advances 
nsin£  the  “IncTime”  procedure.  This  procedure  increments  the  cum  time  cell  on  the  “sim” 
sheet  by  one.  A  subprocedure,  “IncTimeSatChan”  is  also  run  to  advance  time  on  the  “sat 
chan”  sheet.  Fmally,  a  subprocedure,  “UpdateCumSat”  is  run  to  update  the  cumulative 
number  of  satellites  being  used  on  the  “sat  chan”  sheet.  This  procedure  adds  the  satellites 
in  the  sat  chan  column  to  the  cumulative  number  of  channels  used  and  puts  the  new  value 
in  the  cum  sat  column. 
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0.  Start 


1 .  Reset 


2.  To  step  12 
Generate  New  Msg 


Figure  51.  Main  Program  Flow  Chart)  Steps  1  — 11. 
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5.  Execute  the  Decrement  Time  Process 

Step  5  initiates  the  decrement  time  process  at  Step  21 .  After  execution  of  this 
procedure  continue  with  Step  6. 


6.  Stop  Transmission 

Step  6  determines  whether  the  maximum  time  for  each  message  in  the 
retransmission  queue  has  been  reached.  The  procedure  compares  the  time  to  stop  trying 
column  with  the  current  time.  If  the  times  are  equal,  this  procedure  will  stop  the  message 
process.  The  following  is  an  extract  from  the  procedure. 

time  =  Sheets("sim").Cells(9, 2).  Value 
Slieets("Retrans").Select 
maxRow  =  Cells(2, 2).  Value 
For  chkRow  =  srow  To  maxRow 

If  (time  =  Cells(chkRow,  mxTmCol).  Value)  And  (Cells(chlcRow,  stcol).  Value  =  "w") 
Then 

Cells(chkRow,  stcol).  Value  =  "stop" 
origMsgNum  =  Cells(chkRow,  2).  Value 
MsgNumRow  origMsgNum,  origMsgRow,  msgtime 
Sheets("sim").  Select 

If  CeUs(origMsgRow,  34). Value  =  "si"  Then 

Cells(origMsgRow,  34).Value  =  "stopG"  'this  puts  stop  but  msg  good 
in  status  col  on  "sim" 

'indicates  ack  msg  jammed 

Sheets("Retrans").Cells(chkRow,  stcol).  Value  =  "stopG" 

Else 

Cells(origMsgRow,  34). Value  =  "stopNG"  'puts  stop  msg  never  made 
it 

Sheets("Retrans").Cells(chkRow,  stcol).  Value  =  "stq)NG" 

End  If 

Sheets("Retrans").Select 

End  If 

Next 


7.  Execute  the  Retransmit  Message  Process 

Step  7  initiates  the  retransmit  message  process  at  Step  38.  After  execution  of  this 
procedure,  continue  with  Step  8. 
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8.  Execute  the  New  Message  Process 

Step  8  invokes  the  new  message  process  at  Step  12.  After  the  execution  of  this 
process,  continue  with  Step  9. 


9.  Sort 

Step  9  sorts  the  retransmission  queue  by  the  priority  (type)  of  the  message  and 
then  by  time.  This  is  done  to  ensure  that  the  higher  priority  messages  get  retransmitted 
first.  The  following  is  an  extract  of  this  procedure. 


Sheets("Retrans").Select 
startrow  =  7 

sortrow  =  Cells(2, 2).  Value 
If  sortrow  <  100  Then 

Range("A7:V107").Select 

Selection.sort  Keyl:=Range("C7"),  Order  1  :=xlAscending,  Key2:=Range 
("A7"),  Order2:=xIAscending,  Key3:=Range("B7"),  Order3:=  _ 
xlAscending,  Header:=xlGuess,  OrderCustom"!,  MatchCase:=False  _ 

,  Orientation:=xlTopToBottom 
Elself  sortrow  <  200  Then 

Range("A7:V207").Select 

Selection.sortKeyl:=Range("C7"),  Orderl:=xlAscending,  Key2:=Range 
("A7"),  Order2:=xlAscending,  Key3:=Range("B7"),  Qrder3:=  _ 
xlAscending,  Header:=xlGuess,  OrderCustom:=l,  MatchCase:=False  _ 

,  Orientation:=xlTopToBottom 
Elself  sortrow  <  300  Then 

Range("A7:V307").Select 


10.  Next  Iteration 

Step  10  checks  to  see  if  there  is  another  time  step  increment.  If  there  is,  go  to 
Step  3;  otherwise,  continue  to  Step  11. 


11.  End 

This  step  ends  the  program. 
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12.  New  Message  Process 

Step  12  begins  the  new  message  process.  The  first  procedure  inside  this  process  is 
the  new  message  generation  procedure,  “Msg_Input_Lat_Long_Type”.  This  procedure 
generates  new  messages  for  the  current  time  step. 

13.  Number  of  Message  Types 

Step  13  starts  the  “For”  loop  that  occurs  for  each  message  type  in  the  generate 
new  message  procedure. 

14.  Generate  Number  of  Messages  in  Message  Type 

For  each  message  type  a  Poisson  random  number  is  generated  to  determine  the 
number  of  messages  that  will  occur  during  that  time  step.  The  Poisson  random  nmnber 
generator  uses  the  lambda  values  that  were  input  by  the  user.  These  lambda  values  are  the 
average  rate  of  occurrence  for  the  message  type.  The  following  is  an  extract  of  the 
Poisson  random  number  generator  procedmre: 


theNum  =  0 
a  =  E3q)(-lambda) 
b=  1 

Do  Until  Done 

UnifRand  UnifRnd,  stream 
b  =  b  *  UnifRnd 
If  b  <  a  Then 

Done  =  True 

Else 

theNum  =  theNum  +  1 

End  If 

Loop 

mdNum  =  theNum 


15.  Number  of  Messages  for  each  Message  Type 

A  ‘Tor”  loop  is  executed  for  each  message  of  every  message  type. 
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Figure  52.  New  Message  Flow  Chart,  Steps  12  —  20. 
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16.  Generate  New  Message  Information 

Several  Uniform  random  numbers  are  generated  to  create  the  following  data  for 
the  new  message: 


•  a  new  location  for  the  sender  and  receiver;  the  sender  and  receiver  type, 
which  indicates  the  type  of  equipment  the  sender  and  receiver  use  to 
counteract  jamming; 

•  the  number  of  bits,  which  corresponds  to  the  appropriate  type  of  message; 

•  to  determine  if  the  message  is  transmitted  using  the  low  or  medium  data 
rate  if  the  message  is  transmitted  using  satellites  “3”  or  “4”. 

The  following  is  an  extract  of  the  new  message  generation  procedure. 


Cells(mum,  tmcol).  Value  =  Cells(9, 2).  Value 
NewMsg  mum 

Cells(mum,  mtcol).Value  =  msgType 
UnifRand  mdNum,  stream 

Cells(mum,  latcol).Value  =  Int((upboimdLat  -  lowboundLat  +  1)  *  mdNum  + 
lowboundLat) 

UnifRand  mdNum,  stream 

Cells(mum,  latcol  +  1). Value  =  Int((upbndLong  -  lowbndLong  +  1)  *  mdNum  + 
lowbndLong) 

UnifRand  mdNum,  stream 

Cells(mum,  latcol  +  2). Value  =  Int((upboundLat  -  lowboundLat  +  1)  *  mdNum  + 
lowboundLat) 

UnifRand  mdNum,  stream 

Cells(mum,  latcol  +  3). Value  =  Int((upbndLong  -  lowbndLong  +  1)  *  mdNum  + 

'  lowbndLong) 

UnifRand  mdNum,  stream 

Cells(mum,  tcol).Value  =  Int((upbndStype  -  lowbndStype  +  1)  *  mdNum  + 
lowhndStype) 

UnifRand  mdNum,  stream 

Cells(mum,  tcol  +  1).  Value  =  Int((upbndRtype  -  lowbndRtype  +  1)  *  mdNum  + 
lowlmdRtype) 

Cells(mum,  tcol  +  2).  Value  =  msgType  'priority  is  same  as  msgtype  for  now 
UnifRand  mdNum,  stream 

Cells(mum,  mcol  +  2). Value  =  Int((upbits  -  lowbits  +  1)  *  mdNum  +  lowbits) 
UnifRand  mcfi^um,  stream 

Cells(mum,  12). Value  =  mdNum  'puts  random  #  in  "sim"  active  row  rand  #  col  (12) 
ack  =  False  'indicates  this  is  not  an  auto  acknowledgment  message 
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17.  Execute  Process  Information 


Step  17  initiates  the  “process  info”  procedure  at  Step  46,  the  final  operation  for 
processing  messages.  After  execution  of  this  process,  continue  with  Step  18. 

18.  Next  Message 

Step  18  iterates  to  the  next  message  for  the  specified  message  type  firom  Step  15. 
Once  the  required  number  of  messages  for  the  designated  message  type  have  been 
processed,  go  to  Step  19. 

19.  Next  Message  Type 

Step  19  returns  the  program  to  Step  13  if  there  is  another  message  type  to  process, 
otherwise  continue  with  Step  20. 


20.  End  New  Message 

Step  20  ends  the  generate  new  message  process.  Program  control  is  returned  to 
either  Step  2  or  Step  8. 

21.  Decrement  Time 

Step  21  begins  the  Decrement  Time  process  which  decrements  time  on  the  “time” 
sheet.  This  is  done  to  monitor  the  time  remaining  to  transmit  a  message.  The  following  is 
an  excerpt  fi-om  the  “DecrTime”  procedure,  which  encon5)asses  Steps  21  -  37. 


Fot  satNiim  =  1  To  maxSat 

For  chan  =  1  To  maxChan 

If  Cells(crow,  tmcol).  Value  o  ""  Then 

timeLeft  =  Cells(crow,  tmcol).Value 
If  timeLeft  -  1  <=  0  And  (Cells(crow,  7  +  (satNum  -  1)  *  4)  o 
curTime)  Then 

'clears  the  contents  of  the  channel  W  has  completed  its  transmission 
'a  bl^  row  means  it  is  free  to  transmit 
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'the  2nd  part  ensures  that  an  ack  msg  added  does  not  deer  time  in  this  time  step 
msgNum  =  Cells(crow,  mcol  *  satNum). Value 
MsgNumRow  msgNum,  MsgRow,  msgtime 
DecCumSat  satNum 
UpdateCumSat  prevRow 
Sheets("time").Select 

Set  rl  =  Range(Cells(crow,  4  +  (satNum  - 1)  *  4),  Cells(crow, 
5  +  (satNum  -  1)  *  4)) 

Set  r2  =  Range(Cells(crow,  6  +  (satNum  - 1)  *  4),  Cells(crow, 
7  +  (satNum  - 1)  *  4)) 

Set  theRng  =  Union(rl,  r2) 
theRng.Select 
Selection.ClearContents 
Sheets("sim").Select 

If  CeIls(MsgRow,  upCol).Value  =  satNum  Then 
MsgSent  msgNum,  curRow,  stream 

End  If 

Sheets("time").SeIect 

Else 

Cells(crow,  tmcol). Value  =  timeLeft  -  1 

End  If 

End  If 

crow  =  crow  +  1 

Next 

crow  =  ink 
tmcol  =  tmcol  +  4 

Next 


22.  Decrement  Time  for  all  Satellites 

Step  22  starts  the  “For”  loop  that  decrements  time  for  each  satelhte. 


23.  Decrement  Time  for  all  Channels 

Step  23  starts  the  inner  ‘Tor”  loop  of  the  “DecrTime”  procedtire.  This  “For”  loop 
occurs  for  every  channel  on  the  current  satellite. 

24.  Channel  Used 

If  the  current  channel  is  transmitting  a  message,  continue  with  Step  25;  otherwise, 
go  to  Step  35. 
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Figure  53.  Decrement  Time  Fiow  Chart,  Steps  21  -  37. 


25.  Time  Up 

If  time  =  0  on  the  “time”  sheet,  continue  with  Step  26;  otherwise,  go  to  Step  34. 
This  initiates  the  process  of  ending  the  transmission  of  the  message  and  begins  the 
transmission  of  an  acknowledgment  message  if  one  is  necessary. 

26.  Decrement  Cumulative  Number  of  Satellites 

Step  26  decrements  the  number  of  users  on  the  satellite  channel  (“sat  chan”  sheet) 
by  the  procedure  “DecCumSat”.  The  following  is  an  extract  of  this  procedure. 


Sheets("sat  chan").Select 

Done  =  False 

tmcol  =  1 

num=  1 

cnum  =  2 

curRow  =  Cells(9, 2).  Value 
Do  Until  Done 

If  num  =  SatNumToDec  Then 

dec  =  Cells(curRow,  cnum  +  (num  - 1)  *  2). Value 
Cells(curRow,  cnum  +  (num  - 1)  *  2).Value=dec  - 1 
Done  =  True 

Else 

num  =  num  +  1 

End  If 

Loop 


After  execution  of  this  procedure,  the  message  is  also  removed  firom  the  “time” 
sheet  as  indicated  in  the  extract  of  the  “DecrTime”  procedure  in  Step  21 . 


27.  IfUpUnkSateUite 

If  the  current  channel  is  on  the  uplink  satellite,  continue  with  Step  28;  otherwise, 
go  to  Step  35.  This  is  to  ensure  that  Steps  28  -  33  are  executed  only  for  the  uplink 
portion  of  the  transmission. 
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28.  Message  Sent 

The  procedure,  “MsgSent”  initiates  the  process  to  annotate  that  a  message  has 
been  sent.  The  following  excerpt  from  this  procedure  includes  Steps  29-33. 


If  Cells(MsgRow,  AckCol). Value  o  "ack"  Then  If  Cells(MsgRow,  AckCol).Value  o  "ack" 
Then 

SendAck  MsgRow,  curRow,  stream 
GoodMsgNum  msgNum,  goodMsgRow,  status 
Sheets("sim").Select 
Cells(curRow,  l).Value  =  time 
If  status  =  "good"  Then 

Cells(MsgRow,  cnum). Value  =  "si"  'indicates  msg  went  thru  "good"  need  to 
wait  to  see  if  ack  does 

Else 

Cells(MsgRow,  cnum).  Value  =  "s3"  'indicates  msg  was  partially  jammed 

End  If 

Else 

Cells(MsgRow,  cnum).Value  =  "s2"  'indicates  ack  went  thru 
origMsgNum  =  msgNum  *  (-1) 

MsgNumRow  origMsgNum,  origMsgRow,  msgtime 
GoodMsgNum  origMsgNum,  goodMsgRow,  status 
If  status  =  "good"  Then 

'this  will  take  the  original  message  of  the  Retrans  queue 
'othawise  the  msg  will  stay  on  the  retrans  queue  to  be  retransmitted 
'which  is  what  should  happen  if  the  message  was  jammed  at  all 
Sheets("sim").Cells(origMsgRow,  cnum).Value  =  "s" 

TakeOffRetrans  origMsgNum 

End  If 

End  If 


29.  Acknowledgment  Message 

If  the  message  is  an  acknowledgment  message.  Step  29  sends  control  to  Steps  30  - 
31;  otherwise,  go  to  Steps  32  -  33. 


30.  Original  Message  “good” 

If  the  original  message  was  transmitted  “good”,  continue  with  Step  31;  otherwise, 
go  to  Step  35. 
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31.  Remove  Message 

Step  31  removes  the  message  from  the  retransmission  queue.  The  following 
excerpt  is  from  the  “TakeOfiERetrans”  procedure: 


Sheets("Retrans").Select 
Do  Until  Done 

If  msgNum  =  Cells(mum,  mcol). Value  Then 
Cells(minn,  stcol).Value  =  "good" 
Done  =  True 

Else 

mum  =  mum  +  1 

End  If 
Loop 


32.  Send  Acknowledgment 

If  the  message  was  not  an  acknowledgment  message,  then  execute  the  “SendAck” 
procedure.  This  procedure  copies  the  information  from  the  “time”  sheet  to  the  “sim” 
sheet,  but  the  original  sender  is  now  the  receiver  and  the  receiver  is  now  the  sender.  The 
acknowledgment  message  indicates  whether  the  original  message  was  a  “good”  message 
or  a  partially  jammed,  “pj”  message.  The  following  is  an  excerpt  from  the  “SendAck” 
procedure. 


priOTity  =  4 

AckBits  =  2000  'set  num  of  bits  to  transmit  acknowledgment 
Sheets("sim").Select 

Cells(curRow,  1).  Value  =  Cells(9, 2).  Value 

mtyp  =  Cells(MsgRow,  3).Value 

slat  =  Cells(MsgRow,  5).  Value 

slon  =  Cells(MsgRow,  6).  Value 

rlat  =  Cells(MsgRow,  7). Value 

rlon  =  Cells(M^Row,  8). Value 

styp  =  Cells(MsgRow,  9).  Value 

rtyp  =  Cells(MsgRow,  10).  Value 

RandNum  =  Cells(MsgRow,  12).  Value 

jtyp  =  Cells(MsgRow,  27).Value  'indicates  if  original  msg  was  "g"  or"pj" 

Cells(curRow,  3).  Value  =  mtyp 

Cells(curRow,  5).  Value  =  rlat 

Cells(curRow,  6).  Value  =  rlon 

Cells(curRow,  7).  Value  =  slat 
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Cells(curRow,  8).  Value  =  slon 
Cells(curRow,  9).Value  =  rtyp 
Cells(curRow,  10).  Value  =  styp 
Cells(curRow,  4).Value  =  AckBits 
Cells(curRow,  1 1).  Value  =  priority 
Cells(curRow,  12).  Value  =  RandNum 
Cells(curRow,  28).Value  =  "ack" 
Cells(curRow,  30).  Value =jtyp 
mnum  =  Cells(MsgRow,  2).Value 
Cells(curRow,  2). Value  =  mnum  *  (-1) 
ack  =  True 

processinfo  curRow,  stream,  mtyp,  ack 
NewRowSim  curRow 


33.  Execute  Process  Information 

Step  33  initiates  the  “process  info”  procedure,  the  final  operation  for  processing 
messages,  at  Step  46.  After  execution  of  Steps  46  -  63,  continue  with  Step  34. 

34.  Decrement  Time 

This  step  decrements  the  time  by  one  time  unit  if  the  current  message  on  the 
“time”  sheet  is  still  transmitting. 

35.  Next  Channel 

If  there  is  another  channel  on  the  current  satellite,  go  to  Step  23;  otherwise,  go  to 
Step  36. 

36.  Next  Satellite 

If  there  is  another  satellite  to  be  processed,  return  to  Step  22;  otherwise,  go  to 
Step  37. 

37.  End  Decrement  Time 

Step  37  ends  the  Decrement  Time  process  and  control  is  returned  to  Step  5. 
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38.  Retransmit  Message  Process 

Step  38  initiates  the  Retransmit  Message  process.  This  process  retransmits 
messages  that  are  in  the  retransmission  queue. 

39.  Check  for  Retransmission 

Step  39  begins  the  “RetransMsg”  procedure  with  a  “Do”  loop  for  retransmitting 
messages.  This  “Do”  loop  occurs  for  all  messages  in  the  retransmission  queue.  The 
following  is  an  excerpt  from  this  procedure. 


Sheets("sim").Select 

time  =  Cells(curRow,  tmsim).  Value 

Done  =  False 

Sheets("Retrans").Select 

Do  Until  Done 

retranstime  =  Cells(crow,  tmcol).  Value 

If  (Cells(crow,  tmcol). Value  =  time  And  Cells(crow,  stcol).Value  =  "w")  Or  (Cells(crow, 
tmcol). Value  <  time  And  Cells(crow,  stcol). Value  =  "w")  Then 

the  second  part  of  the  if  statement  indicates  that  there  werent 
available  channels  previously 
AvailChan  AvailCh,  crow,  "Retrans" 

If  AvailCh  Then 

sendRetrans  crow,  curRow,  msgType 
processinfo  curRow,  stream,  msgType,  ack 
Sheets("sim").Select 
If  Cells(curRow,  34).Value  =  "w"  Then 
Sheets("Retrans").Select 
Cells(CTOw,  14).  Value  =  "w" 

Elself  CelIs(curRow,  34).Value  =  "t"  Then 
Sheets("Retrans").Select 
Cells(crow,  14),  Value  =  "t" 

End  If 

NewRow  curRow 
crow  =  CTOW  +  1 

Else 

CTOW  =  crow  +  1 

End  If 

Elself  Cells(CTOw,  tmcol).  Value  =  ""  Then 
Done  =  True 

Else 

crow  =  crow  +  1 

End  If 

Sheets("Retrans").Select 

Loq) 
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Figure  54.  Retransmit  Message  Flow  Chart,  Steps  38  -  46. 


40.  Time  to  Retransmit 


If  the  time  in  the  time  to  transmit  column  is  less  than  or  equal  to  the  current  time, 
and  the  status  column  in  the  “Retrans”  sheet  indicates  waiting,  “w”,  continue  with 
Steps  41  to  43;  otherwise,  go  to  Step  44. 


41.  Available  Channel 

If  there  is  an  available  channel  for  both  the  uplink  and  the  downlink  on  the  current 
message,  continue  with  Steps  42  and  43.  Otherwise,  if  there  is  not  an  available  channel, 
go  to  Step  44.  The  following  is  an  extract  of  the  “AvailChan”  procedure. 


If  sheetnm  =  "sim"  Then 
cnumUp  =14 

Elself  sheetnm  =  "Retrans"  Then 
cnumUp=  16 

End  If 

cnumDown  =  17 
Sheets(sheetnni).Select 

satNumUp  =  Cells(curRow,  cnnmUp).  Value  'identifies  the  current  sat  uplink 

satNumDown  =  Cells(curRow,  cnumDown).Value 

Sheets("sat  chan").Select 

satChRow  =  Cells(9, 2).  Value 

satMax  =  Cells(6, 3).  Value 

Do  Until  Dcme 

If  Cells(ChRow,  cnum).  Value  =  satNumUp  Then 
satMaxCh  =  Cells(satRow,  cnum).  Value 
If  Cells(satChRow,  cnum).  Value  <  satMaxCh  Then 
AvailCh  =  True 
Done  =  True 

Else 

AvailCh  =  False 
Done  =  True 

End  If 

Else 

cnum  =  cnum  +  2 

End  If 

Loop 

If  satNumUp  o  satNumDown  And  AvailCh  =  True  Then 
Done  =  False 
cnum  =  3 
Do  Until  Done 

If  Cells(ChRow,  mum).  Value  =  satNumDown  Thm 
satMaxCh  =  Cells(satRow,  cnum).  Value 
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Loop 

End  If 


If  Cells(satChRow,  cnum).  Value  <  satMaxCh  Then 
AvailCh  =  True 
Done  =  True 

Else 

AvailCh  =  False 
Done  =  True 

End  If 

Else 

cnum  =  cnum  +  2 

End  If 


42.  Send  Retransmission 

Step  42  attempts  to  resend  the  message  by  calling  the  “SendRetrans”  procedure. 
Information  from  the  retransmission  queue  is  copied  into  the  “sim”  sheet.  The  following  is 
an  extract  from  the  procedure. 


Sheets("Retrans").Select 

Set  rl  =  Range(Cells(RetRow,  2),  Cells(RetRow,  3)) 

Set  r2  =  Range(Cells(RetRow,  4),  Cells(RetRow,  5)) 

Set  r3  =  Range(Cells(RetRow,  6),  Cells(RetRow,  7)) 

Set  r4  =  Range(Cells(RetRow,  8),  Cells^etRow,  9)) 

Set  r5  =  Range(Cells(RetRow,  10),  Cells(RetRow,  1 1)) 

Set  retRng  =  Union(rl,  r2,  r3,  r4,  r5) 

retRng.Select 

Selection.copy 

Sheets("sim").Select 

Set  srl  =  Range(Cells(curRow,  2),  Cells(curRow,  3)) 

Set  sr2  =  Range(Cells(curRow,  4),  Cells(curRow,  5)) 

Set  sr3  =  Range(Cells(curRow,  6),  CeIls(curRow,  7)) 

Set  sr4  =  Range(Cells(curRow,  8),  Cells(curRow,  9)) 

Set  sr5  =  Range(Cells(curRow,  10),  Cells(curRow,  1 1)) 

Set  simRng  =  Union(srl,  sr2,  sr3,  sr4,  sr5) 
simRng.Select 

Selection.PasteSpecial  Paste:=xlValues,  Operation  :=xlNone, 

SkipBlanks:=False,  Transpose:=False 

time  =  Cells(9, 2).  Value 

Sheets("Retrans").Select 

RandNum  =  Cells(R^ow,  18).  Value 

Sheets(“sim”).Select 

Cells(curRow,  12).  Value  =  RandNum 

Sheets("Retrans").Select 

Cells(RetRow,  14).  Value  =  "r" 

Sheets("sim").Select 
Cells(curRow,  29).  Value  =  V 
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Cells(curRow,  1).  Value  =  time 

mtype  =  Cells(curRow,  3).  Value  'msg  type;  value  needed  for  processing  info 

43.  Execute  Process  Information 


Step  43  initiates  the  “process  info”  procedure,  the  final  operation  for  processing 
messages,  at  Step  46.  After  execution  of  Steps  46  -  63  continue  with  Step  44. 


44.  Loop 

Step  44  determines  if  there  is  another  message  in  the  retransmission  queue.  If 
there  is  another  message  in  the  queue,  go  to  Step  39;  otherwise,  continue  to  Step  45. 


45.  End  Retransmission 

Step  45  ends  the  retransmission  procedure  and  control  is  returned  to  Step  7. 


46.  Process  Information 

Step  46  begins  the  Process  Information  procedure. 


47.  Satellite  Type 

Step  47  identifies  the  satellites  that  will  transmit  the  messages.  The  procedure 
con^ares  the  location  of  the  sender  and  receiver  to  the  appropriate  satellite  that  will 
transmit  the  message.  The  following  is  an  extract  fi-om  the  satellite  type  procedure. 

maxSat  =  Sheets("sat  chan").Cells(2, 4).Value 
Sheets("sim").Select 

If  ((Cells(curRow,  cstrt).  Value  >=  0  And  Cells(curRow,  cstrt).  Value  <  90)  Or  Cells(curRow, 
cstrt).  Value  =  360)  Then 

If  Sheets("sat  chan").CelIs(rstat,  cstat).  Value  =  True  Then 
Sheets("sim").CeIIs(curRow,  satUp).Value  =  1 
Elself  maxSat  >  4  Then 

If  Sheets("sat  chan").CeIIs(rstat,  cstat  +  8). Value  =  True  Thai 
Sheets("sim").Cells(curRow,  satUp).Value  =  5 

End  If 

Else 

Sheets("shn").Cells(curRow,  satUp).  Value  =  1 

End  If 
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Elself  (Cells(curRow,  cstrt).  Value  >=  90  And  Cells(curRow,  cstrt).  Value  <  180)  Then 
If  SheetsC'sat  chan").Cells(rstat,  cstat  +  2). Value  =  True  Then 
Sheets("sim").Cells(curRow,  satlJp).Value  =  2 
Elself  maxSat  >  5  Then 

If  SheetsC'sat  chan").Cells(rstat,  cstat  +  10).Value  =  True  Then 
Sheets("sini").Cells(curRow,  satUp).Value  =  6 

End  If 

Else 

Sheets("sini").Cells(curRow,  satlJp).Value  =  2 

End  If 

Elself  (Cells(curRow,  cstrt).  Value  >=  180  And  Cells(curRow,  cstrt).Value  <  270)  Then 
If  SheetsC'sat  chan").Cells(rstat,  cstat  +  4).  Value  =  True  Then 
Sheets("sini").Cells(curRow,  satUp).  Value  =  3 
Elself  maxSat  >  6  Then 

If  SheetsC'sat  chan").Cells(rstat,  cstat  +  12).Value  =  True  Then 
Sheets("sim").Cells(curRow,  satUp).Value  =  7 

End  If 

Else 

Sheets("sini").Cells(curRow,  satUp).  Value  =  3 

End  If 

Elself  (Cells(curRow,  cstrt).Value  >=  270  And  Cells(curRow,  cstrt).Value  <  360)  Then 
If  SheetsC'sat  chan").Cells(rstat,  cstat  +  6).  Value  =  True  Then 
Sheets("sim").Cells(curRow,  satUp). Value  =  4 
Elself  maxSat  >  7  Then 

If  SheetsC'sat  chan").Cells(rstat,  cstat  +  14).  Value  =  True  Then 
Sheets("sim").Cells(curRow,  satUp). Value  =  8 

End  If 

Else 

Sheets("sim").Cells(curRow,  satUp).  Value  =  4 

End  If 

End  If 

If  ((CelIs(curRow,  cstp). Value  >=  0  And  Cells(curRow,  cslp).Value  <  90)  Or  Cells(curRow, 
cstp).  Value  =  360)  Then 

If  SheetsC'sat  chan'').Cells(rstat,  cstat).  Value  =  True  Then 
Sheets("sim").Cells(curRow,  satDown).Value  =  1 
Elself  maxSat  >  4  Then 

If  SheetsC'sat  chan").Cells(rstat,  cstat  +  8).  Value  =  True  Then 
Sheets("sim").Cells(curRow,  satDown).  Value  =  5 

End  If 

Else 

Sheets("sim").Cells(curRow,  satDown).Value  =  1 

End  If 


48.  Copy  Informatioii 

This  step  copies  the  eqtiation  information  from  the  previous  row  to  the  current 
row.  The  columns  copied  are  data  rate,  uplink  and  crosslink  (if  ^plicable),  transmission 
rate,  and  time  to  transmit. 
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Figure  55.  Process  Information  Flow  Chart,  Steps  47  -  51. 
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49.  Location  and  Type 

Step  49  inputs  the  sender  and  receiver  location  and  type  of  user  from  the  “sim” 
sheet  to  the  “type”  and  “jammer”  sheets. 


50.  Percentage  Jammed 

The  procedure  “perc”  is  implemented  in  Step  50.  This  procedure  inputs  the  anti¬ 
jam  capability  proportion  into  the  “jammer”  sheet.  The  procedure  goes  to  the  “type”  sheet 
and  matches  the  type  of  user  to  the  anti-jamming  capability  and  records  this  information 
on  the  “jammer”  sheet.  The  “jammer”  sheet,  previously  discussed,  provides  additional 
information  on  how  this  value  is  used.  The  following  is  an  extract  from  this  procedure. 


ctypenum  =  9 
Sheets(”sim”).Select 

satJamType  =  Cells(curRow,  ctypenmn).Value 
Sheets(”type").Select 

rtypeup  =  3  'row  for  type  sender  from  "type”  sheet 
rtypedown  =  5  'row  for  type  receiver  from  "type" 
ctypeinit  =  4 

Rtype  =  9  'first  row  to  compare  type  # 

ctype  =  1  *col  #  to  compare  type  # 

cper  =  2  'percentage  col  to  return  %  value 

jrow  =  9  'from  "jammer"  sheet  jam  row 

jcolup  =  7  'from  "jammer"  sheet  jam  col  for  up 

jcoldown  =  16  'from  "jammer"  sheet  jam  col  from  down 

Done  =  False 

Do  Until  Done 

If  Cells(Rtype,  ctype). Value  =  Cells(rtypeup,  ctypeinit).Value  Then 
theperc  =1-  Cells(Rtype,  cper).  Value 
Sheets("jammer").Select 
Cells(jrow,  jcolup).  Value  =  theperc 
UnifRand  mdNum,  stream 

CellsQrow,  12).Value  =  mdNum  'puts  random  #  in  "jammer"  rand  #  col 
Done  =  Tme 

Else 

Rtype  =  Rtype  +  1 

End  If 

Loop 

Sheks("type").Select 

If  Cells(rtypeup5  ctypeinit). Value  =  Cel^rtypedown,  ctypeinit).  Value  Then 
Sheets("jammer").Select 

Cells(jrow,  jcoldown).  Value  =  Cells(jrow Jcolup).  Value 
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Else 


Done  =  False 
Rtype  =  9 
Do  Until  Done 

If  Cells(Rtype,  ctype).  Value  =  Cells(rtypedown,  ctypeinit).  Value  Then 
theperc  =  1-  Cells(Rtype,  cper).  Value 
Sheets("jammer").Select 
Cells(jrow,  jcoldown).  Value  =  theperc 
UnifRand  mdNum,  stream 

Cells(jrow,  21). Value  =  mdNum  'puts  random  #  in  "jammer"  rand  # 
col 

Done  =  True 

Else 

Rtype  =  Rtype  +  1 

End  If 

Loop 

End  If 


51.  Available  Channel 

This  step  determines  whether  there  is  an  available  channel  on  the  satellite(s)  for 
both  the  uplink  and  downlink  (if  necessary).  It  selects  the  available  channel  number  from 
the  “time”  sheet  and  returns  the  channel  number,  the  row  for  the  channel  number,  and  the 
column  number  for  the  available  satelhte.  It  utilizes  this  chaimel  next  on  the  satellite  that 
is  transmitting.  The  following  is  an  extract  from  this  procedure. 


Sheets("time").Select 
Do  Until  SatDone 

If  satUp  =  CeIIs(satRow,  cnum).  Value  Then 
Done  =  False 
Do  Until  Done 

If  Cells(crow,  cnum). Value  =  ""  Then 

ChNumUp  ==  Cells(crow,  ChCol).  Value 
RChNumUp  =  crow 
ColChUp  =  cnum 
Done  =  Tme 

Else 

crow  =  crow  +  1 

End  If 

Loop 

SatDone  =  True 

Else 

cnum  =  oium  +  4 

End  If 

Loq) 
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If  satUp  o  satDown  Then 
SatDone  =  False 
cnum  =  4 
crow=  12 
TimeCol  =  5 
Do  Until  SatDone 

If  satDown  =  Cells(satRow,  cnum).  Value  Then 

Done  =  False 

Do  Until  Done 

If  Cells(crow,  cnum).  Value  =  ""  Then 

ChNumDown  =  Cells(crow,  ChCol). Value 
RChNumDown  =  crow 
ColChDown  =  cnum 
Done  =  True 

Else 

crow  =  crow  +  1 

End  If 

Loop 

SatDone  =  True 

Else 

cnum  =  cnum  +  4 

End  If 


Loop 

Else 

ChNumDown  =  ChNumUp 
RChNumDown  =  RChNumUp 
ColChDown  =  ColChUp 

End  If 


52.  Check  Available  Channel 

If  there  is  an  available  channel,  go  to  Step  53;  otherwise,  go  to  Step  56.  Step  41 
provides  an  extract  of  the  “AvailChan”  procedure. 

53.  Communication  Occurred 

Step  53  calls  the  “CommOccur”  procedure.  This  procedure  determines  whether 
the  transmission  was  “good”,  meaning  that  no  jamming  occurred;  “pj”  indicating  that  the 
message  was  partMy  jammed;  or  “j”  indicating  that  the  message  was  con^letely  jammed. 
The  following  is  an  excerpt  from  this  procedure. 

Sheets("jammer").Select 

If  Cells(crow,  ccoldown).  Value  =  2  Then  *2  indicates  jam  capability  was  100% 

JamTotal  curRow,  msgType,  ack 

Elself  Cells(crow,  ccoldown).  Value  =  1  Then  T  indicates  jam  c^  was  partial 
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JamPartial  curRow,  msgType,  ack 

Else 

JamNone  curRow,  msgType,  ack 

End  If 


54.  Check  for  Jam  Total 

Step  54  checks  to  see  if  the  message  was  totally  jammed.  If  it  was,  continue  with 
Step  55;  otherwise,  go  to  Step  59. 


55.  Jam  Total 

Step  55  puts  a  “j”  on  the  “sim”  sheet  to  indicate  that  the  message  was  completely 
jammed.  The  following  is  an  extract  from  the  “JamTotal”  procedure. 

Sheets("sim").Select 

Cells(curRow,jcol). Value  =  "j"  'puts  on  "sim"  sheet  wiiat  jamming  occurred  to  msg;  value  will 
never  change 

Cells(curRow,  ccol).  Value  =  "j"  'indicates  the  satellite  is  jammed 
timLastTry  =  CelIs(curRow,  1).  Value 
If  (Not  ack)  Then 

WaitTimePri  \raitTime,  msgType,  maxWaitTime 

NeedToRetrans  waitTime,  curRow,  maxWaitTime,  timLastTry  'puts  message  in 
retransmission  queue  wait 

End  If 


56.  Wait  Time  by  Priority 

Step  56  calls  the  procedure  “WaitTimePri”  to  determine  the  amount  of  time  the 
message  should  wait  before  trying  to  retransmit.  The  amoxmt  of  time  is  determined  by  the 
type  of  message.  This  procedure  goes  to  the  “type”  sheet  and  matches  the  message  type 
value  to  the  wait  time  and  maximum  wait  time  values.  The  wait  time  value  is  the  amount 
of  time  the  sender  will  wait  for  an  acknowledgment  message  before  trying  to  send  the 
message  again.  The  maximum  wait  time  value  is  the  maximum  amoxmt  of  time  the  sender 
is  willing  to  continue  trying  the  message  without  receiving  an  acknowledgment  message. 
The  following  is  an  extract  of  this  procedure. 
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Avail 

Chan 


yes 


60.  Jam 
None 


_ _ 

55.  Jam  Total 


56.  Wait  Time  by 
Priority  of  Msg 


_ Nk _ 

61.  Put  Msg  on 
Channel 


62.  Update  Cum 
Sat;  Increment  # 
of  Users 


Figure  56.  Process  lufonnatioii  Flow  Chart  continued.  Steps  52  —  63. 


Sheets("type").Select 
Done  =  False 
Do  Until  Done 

If  msgType  =  Cells(TRow,  typCol).Value  Then 

waitTime  =  Cells(TRow,  typCol  +  1).  Value 
maxWaitTime  =  Cells(mxTrow,  typCol  +  l).Value 
Done  =  True 

Else 

TRow  =  TRow  +  1 
mxTrow  =  mxTrow  +  1 

End  If 

Loop 

57.  Need  To  Retransmit 


Step  57  calls  the  procedure  “NeedToRetrans”.  This  procedure  puts  the  message 
on  the  retransmission  queue  by  copying  the  message  information  from  “sim”  sheet  into  the 
“Retrans”  sheet.  It  puts  a  “w”  in  the  status  column  in  the  “Retrans”  sheet  to  indicate  that 
the  message  is  waiting  to  be  retransmitted,  or,  if  the  message  was  a  retry  message,  it  puts 
an  “n”  indicating  that  the  message  was  a  retransmission.  The  procedure  also  puts  the 
mnxiTniim  time  to  keep  trying  to  send  the  message  and  the  amount  of  time  to  wait  prior  to 
retransmitting  a  message  into  the  retransmission  queue.  The  following  is  an  excerpt  from 
this  procedure. 


Sheets("Retrans").Select 

RetRow  =  Cells(2, 2).Value 

tmcol  =  12 

mxTmCol  =  20 

Sheets("sim").Select 

curMs^um  =  Cells(curRow,  2).  Value 

MsgNumRow  curMsgNum,  origMsgRow,  origMsgTime  'determines  original  message  row  on 
"sim" 

Sheets("sim").Select 
Cells(curRow,  l).Select 
Selection.copy 
Sheets("Retrans").Select 
Cells(RetRow,  l).Select 

Selection.PasteSpecial  Paste:=xlValues,  Operation:=xlNone,  _ 

SkipBlanks:=False,  Transpose:=False 
Sheets("sim").Select 

Set  srl  =  Range(Cells(curRow,  2),  Cells(curRow,  3)) 
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Set  sr2  =  Range(Cells(ciirRow,  4),  Cells(curRow,  5)) 

Set  sr3  =  Range(Cells(curRow,  6),  CeIls(curRow,  7)) 

Set  sr4  =  Range(Cells(curRow,  8),  Cells(curRow,  9)) 

Set  sr5  =  Range(Cells(curRow,  10),  Cells(curRow,  11)) 

Set  simRng  =  Union(srl,  sr2,  sr3,  sr4,  sr5) 

simRng.Select 

Selection.copy 

Sheets("Retrans").Select 

Set  rl  =  Range(Cells(RetRow,  2),  Cells(RetRow,  3)) 

Set  r2  =  Range(Cells(RetRow,  4),  Cells(RetRow,  5)) 

Set  r3  =  Range(CeIls(RetRow,  6),  CelIs(RetRow,  7)) 

Set  r4  =  Range(Cells(RetRow,  8),  Cells(RetRow,  9)) 

SetrS  =  Range(CeIls(RetRow,  10),  Cells(RetRow,  11)) 

Set  retRng  =  Union(r  1 ,  r2,  r3,  r4,  r5) 
retRng.Select 

Selection.PasteSpecial  Paste:=xlValues,  OperatiQn:=xlNone,  _ 

SkipBlanks:=False,  Transpose:=False 
Cells(RetRow,  tmcol).Value  =  waitTime 
Cells(RetRow,  mxTmCol).Value  =  maxWaitXime 
Sheets("sim").Select 
If  Cells(curRow,  rcol).Value  =  "r"  Then 
status  =  "n" 

RetMsgNumRow  curMsgNum,  RetOrigRow 
Sheets("Retrans").Select 

Cells(RetOrigRow,  15). Value  =  timLastTry  'puts  time  last  tried  to  trans  in  orig  msg  row 
Cells(RetOrigRow,  13).Value  =  Cells(RetOrigRow,  tmcol).Value  +  Cells(RetOrigRow, 
15).  Value 

Else 

status  =  "w" 

Sheets("Retrans").Select 
Cells(RetRow,  15).  Value  =  timLastTry 

Cells(RetRow,  13).  Value  =  Cells(RetRow,  tmcol).Value  +  Cells(RetRow,  15).Value 
Cells(RetRow,  21). Value  =  Cells(RetRow,  l).V^ue  +  Cells(RetRow,  mxTmCol). Value 

End  If 

Sheets("sim").SeIect 
Cells(curRow,  cnum).Value  =  "w" 
satUp  =  CelIs(curRow,  14).Value 
satDown  =  Cells(curRow,  17).Value 
RandNum  =  Cells(curRow,  12).Value 
jam  =  CelIs(curRow,  27).Value 
Sheets("Retrans").Select 
Cells(RetRow,  cRetstat).Value  =  status 
Cells(RetRow,  16).  Value  =  satUp 
Cells(RetRow,  17).  Value  =  satDown 
Cells(RetRow,  18).  Value  =  RandNum 
Cells(RetRow,  1 9).  Value  =  jam 
Cells(2, 2).  Value  =  RetRow  +  1 

58.  Check  for  Jam  Partial 

Step  58  checks  to  see  if  the  message  is  partially  jammed.  If  it  was,  then  continue 
with  Step  59;  otherwise,  go  to  Step  60. 
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59.  Jam  Partial 


Step  59  puts  a  “pj”  on  the  “sim”  sheet  to  indicate  that  the  message  was  partially 
jammed.  Following  this  step,  go  to  Step  61.  The  following  is  an  excerpt  from  this 
procedure. 


Sheets("sim").Select 
Cells(curRow,jcoI).  Value  =  "pj" 

Cells(curRow,  ccol).VaIue  =  "pj"  'indicates  message  is  being  partially  jammed 
timLastTry  =  Cells(curRow,  1).  Value 

PutTransOnChan  curRow  'puts  the  transmission  on  the  channel  that  is  trans  on  "time" 
IncNumUsers  curRow  'increments  #  of  cum  users  on  "sat  chan" 

If  (Not  ack)  Then 

WaitTimePri  waitTime,  msgType,  maxWaitTime 

NeedToRetrans  waitTime,  curRow,  maxWaitTime,  timLastTry  'puts  message  into  the 
retransmission  queue 


End  If 


60.  Jam  None 

If  the  message  was  not  jammed  either  partially  or  completely,  a  “good”  is  placed  in 
the  status  column  on  the  “sim”  sheet.  The  following  is  an  excerpt  from  this  procedure. 


Sheets("sim").SeIect 
CeIls(curRow,jcol).  Value  =  "good" 

Cells(curRow,  ccol). Value  =  "t"  'indicates  message  transmitting  good  on  sim  sheet 
timLastTry  =  Cells(curRow,  1).  Value 

PutTransOnChan  curRow  'puts  the  transmission  on  the  channel  that  is  transmitting  on  sheet 
time 

IncNumUsers  curRow  'ino'ements  the  cumulative  number  of  users  on  a  satellite  on  sheet  sat 
chan 

If  (Not  ack)  Then 

WaitTimePri  waitTime,  msgType,  maxWaitTime 
NeedToRetrans  waitTime,  curRow,  maxWaitTime,  timLastTry 

End  If 


61.  Put  Transmission  on  Channel 

Step  61  puts  the  message  on  the  “time”  sheet  indicating  that  the  message  is 
transmitting  by  the  ‘TutTransOnChan”  procedure.  It  puts  the  message  on  both  the  uplink 
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and  downlink  satellites.  The  procedure  copies  the  message  information  from  “sim”  sheet 
to  the  “time”  sheet. 


AvailChanNum  ChUp,  ChDown,  RChUp,  RChDown,  ColChUp,  ColChDown,  curRow 

Sheets("sitn’').Select 

satUp  =  Cells(curRow,  cnumUp).Value 

satDown  =  Cells(curRow,  cnumDown).Value 

Cells(curRow,  cnumMsg).  Select 

Selection.copy 

Sheets("time").Select 

Cells(RChUp,  ColChUp).Select 

Selection.PasteSpecial  Paste:=xlValues,  Operation:=xlNone,  _ 

SkipBlanks:=False,  Transpose:=False 

Sheets("sim").Select 

curTime  =  CeIls(curRow,  1).  Value 

Cells(curRow,  cnumtm).  Select 

Selection.copy 

Sheets("time").Select 

Cells(RChUp,  ColChUp  +  l).Select 

Selection.PasteSpecial  Paste:=xlValues,  Operation:=xlNone,  _ 

SkipBlanks:=False,  Transpose:=False 

Cells(RChUp,  ColChUp  +  2). Value  =  Cells(RChUp,  ColChUp  +  1). Value 
Cells(RChUp,  ColChUp  +  3).  Value  =  curTime 
If  satUp  o  satDown  Then 

Sheets("sim").Select 

Cells(curRow,  cnumMsg).Select 

Selection.copy 

Sheets("time").Select 

Cells(RChDown,  ColChDown).Select 

Selection.PasteSpecial  Paste;=xlValues,  Operati(Mi:=xlNone,  _ 

SkipBlanks:=False,  Transpose:=False 
Sheets(''sim").Select 
Cells(curRow,  cnumtm).Select 
Selection.copy 
Sheets("time").Select 

Cells(RChDown,  ColChDown  +  l).Select  - 
Selection.PasteSpecial  Paste:=xlValues,  Operation:=xlN<Mie,  _ 

SkipBlanks:=False,  Transpose:=False 

Cells(RChDown,  ColChDown  +  2).Value  =  Cells(RChDown,  ColChDown). Value 
Cells(RChDown,  ColChDown  +  3). Value  =  curTime 

End  If 


62.  Increment  Number  of  Users 

Step  62  updates  the  “sat  chan”  sheet.  The  procedure,  ‘TncNumUsers”,  increments 
the  number  of  users  of  the  satellite  that  is  transmitting  by  one.  It  does  this  for  both  the  up 
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and  downlink  satellites.  After  executing  this  procedure  the  program  goes  to  Step  55.  The 
following  is  an  extract  fi-om  this  procedure. 


SheetsC’sat  chan").Select 
curRow  =  Cells(9, 2).  Value 

AvailChanNum  ChUp,  ChDown,  RChUp,  RChDown,  CChUp,  CChDown,  simRow 

Sheets("sim").Select 

satUp  =  Cells(simRow,  cnumUp).  Value 

satDown  =  Cells(simRow,  cnumDown).  Value 

SheetsC'sat  chan").Select 

prevtime  =  Cells(curRow  - 1,  tmcol).  Value 

Cells(curRow,  tmcol).Value  =  prevtime  +  1  'increases  time  value  in  "sat  chan" 

ColChUp  =  satUp*2 

ColChDown  =  satDown  *  2 

incUp  =  Cells(curRow,  ColChUp).Value 

Cells(curRow,  ColChUp).Value  =  incUp  +  1 

If  satUp  o  satDown  Then 

incDown  =  Cells(curRow,  ColChDown). Value 
Cells(curRow,  ColChDown).Value  =  incDown  +  1 

End  If 

nxtRow  =  curRow  +  1 
Cells(9, 2).  Value  =  ipctRow 
UpdateCumSat  curRow 
UpdateCumSat  nxtRow 


63.  End 

Step  63  ends  “processinfo”.  Program  control  is  returned  to  the  appropriate  step 
(i.e.,  Step  17, 33,  or  43). 
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No.  of  messages  No.  of  messages  No.  of  messages 


No.  of  messages  No.  of  messages  No.  of  messages 


No.  of  messages  No.  of  messages 


No.ofchannelsinus.  No.  of  channels  in  use 


No.  of  channels  in  usi  |  No.  of  channels  in  usi  |  |  No.  of  channels  in  use 


The  number  of  channels  in  use  over  time  for  Case  2  is  shown  in  Figures  82  -  86 


Figure  82.  Case  2:  Priority  1. 


0  3  6  9  12  15  18  21  24  27  30  33  36  39  42  45  48  51  54  57  60 
time 


|a  sat1  ■  sat2  □  sat3  n  sat^ 


Figure  83.  Case  2:  Priority  2. 


Figure  84.  Case  2:  Priority  3. 


Figure  85.  Case  2:  Priority  4. 


0  3  8  9  12  16  18  21  24  27  30  33  38  39  42  45  48  51  54  57  60 
time 


|ei  satl  ■  sat2  □  saiiS  □  sat4  ( 


Figure  86.  Case  2:  All  Messages. 
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#  of  channels  in  use  |  #  of  channels  in 


The  number  of  channels  in  use  over  time  for  Case  3  is  shown  in  Figures  87-91 


Figure  87.  Case  3:  Priority  1. 


Figure  88.  Case  3;  Priority  2. 


Figure  89.  Case  3:  Priority  3. 


Figure  90.  Case  3:  Priority  4. 
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30 
25 
20 
15 
10 
5 
0 

0  3  6  9  12  15  18  21  24  27  30  33  36  39  42  45  48  51  54  57  60 
time 

|b  satl  ■  sat2  □  sat3  a  sa^ 


Figure  91.  Case  3:  All  Messages. 


The  number  of  channels  in  use  over  time  for  Case  4  is  shown  in  Figures  92  -  96. 


Figure  94.  Case  4:  Priority  3. 


Figure  95.  Case  4:  Priority  4. 


Figure  96.  Case  4:  All  Messages. 
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APPENDIX  D.  CUMULATIVE  AND  AVERAGE  NUMBER  OF  MESSAGES 

GENERATED 


The  cumulative  number  of  messages  generated  by  priority  and  the  average  number 
of  messages  generated  over  time  is  shown  by  case  in  Figures  97  -  102. 


|>-4^Priority1  --•-^Pfloriiy2  -A-Pr1ority3  -^Prtor»y4  [ 


Figure  97.  Case  2. 


Figure  99.  Case  3. 


Figure  101.  Case  4. 
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APPENDIX  E.  SUMMARY  STATISTICS 


The  following  are  the  summary  statistics  for  runs  of  the  model. 


Casel 

cum  chan  used  time  60 

#  good  trans 

#  stopG 

sat  1 

sat  2 

sat  3 

sat  4 

total 

pi 

p2 

P3 

p4 

total  pi 

p2 

p3  p4 

total 

run  1 

8 

8 

8 

7 

31 

27 

23 

12 

9 

71  13 

20 

5  5 

43 

run  2 

8 

8 

7 

8 

31 

29 

55 

14 

25 

123  16 

17 

4  4 

41 

run  3 

8 

8 

8 

8 

32 

12 

36 

12 

9 

69  19 

28 

5  9 

61 

run  4 

8 

5 

7 

8 

28 

16 

19 

5 

4 

44  15 

29 

6  15 

65 

runs 

8 

8 

7 

8 

31 

24 

28 

9 

13 

74  19 

18 

6  11 

54 

Mean 

8 

7.4 

7.4 

7.8 

30.6 

21.6 

32.2 

10.4 

12 

76.2  16.4 

22.4 

5.2  8.8 

52.8 

Var 

0 

1.8 

0.3 

0.2 

23 

53.3 

202.7 

12,3 

63 

827.7  6.8 

32.3 

0.7  20.2 

113.2 

stdev 

0 

1.3416 

0.5477 

0.4472 

1.S16SS 

7.30068 

14.237 

3.5071 

7.937 

28.77  2.6077 

5.683  8367  4,494 

10.64 

Casel 

#stopNG 

gd/tot#msgs 

Time  to  transmit  gd  msgs 

sat  1 

sat  2 

sat  3 

sat  4 

total 

pi 

p2 

P3 

p4 

total  pi 

p2 

P3  p4 

total 

run  1 

77 

84 

34 

96 

291 

0.20611 

0.1429 

0.1846 

0.0511 

0.1332  5.9259 

9.087 

11  10.11 

833803 

run  2 

50 

67 

26 

86 

229 

0.28431 

0.3056 

0,2154 

0.1344 

0J308  5.5172 

8.673 

13.5  11.84 

9.12195 

run  3 

75 

86 

26 

93 

280 

0.10345 

0.1935 

0.2143 

0.0492 

0.1275  6,1667 

7.694 

10.9  14.56 

8.88406 

run  4 

61 

103 

34 

98 

296 

0.16162 

0.0941 

0.0806 

0.0211 

0.0796  5.25 

8.632 

15.4  17.75 

9 

runs 

71 

97 

20 

78 

266 

0.192 

0.1538 

0.18 

0.0760 

0.1402  5.7917 

7.821 

14.3  16.77 

9.52703 

Mean 

66.8 

87.4 

28 

90.2 

272.4 

0.18950 

0.1780 

0.1750 

0.0664 

0.1423  5.7303 

8.381 

13.0  14.21 

8.97421 

Var 

126.2 

191.3 

36 

67.2 

7213 

0,00436 

0.0063 

0.0030 

0.0018 

0.0030  0.1271 

0.358  4.03  10.40 

0.18524 

stdev 

11.234 

13.831 

6 

8.1976 

26357 

0.06605 

0.0796 

0.0552 

0.0427 

0.0549  0.3566 

0.598  2.01  3.226 

0.43039 

Case2 

cum  chan  used 

time  60 

#  good  trans 

#  StopG 

sat  1 

sat  2 

sat  3 

sat  4 

total 

pi 

p2 

P3 

p4 

total 

pi 

p2 

p3  p4 

total 

run  1 

16 

13 

16 

16 

61 

18 

36 

10 

22 

86 

22 

23 

10  15 

70 

run  2 

16 

16 

8 

12 

52 

41 

51 

15 

44 

151 

19 

15 

7  17 

58 

run  3 

12 

9 

12 

16 

49 

29 

46 

18 

40 

133 

15 

25 

6  11 

57 

run  4 

8 

16 

8 

11 

43 

35 

55 

12 

48 

150 

19 

27 

4  19 

69 

runs 

16 

13 

9 

13 

51 

35 

30 

14 

27 

106 

16 

18 

11  8 

53 

Mean 

13.6 

13.4 

10.6 

13.6 

51.2 

31.6 

43.6 

13.8 

36.2 

1253 

183 

21.6 

7.6  14 

61.4 

Var 

12.8 

8.3 

11.8 

5.3 

42.2 

75.8 

108,3 

9.2 

125.2 

811.7 

7.7 

24.8 

8.3  20 

583 

stdev 

3.5777 

2.8810 

3.4351 

2.3022 

6.49615 

8.70632 

10.407 

3.0332 

11.189 

28.490 

2.7749 

4.980  2.88  4.472 

7.63544 

Case2 

#stopNG 

gd/tod^msgs 

Time  to  transmit  gd  msgs 

sat  1 

sat  2 

sat  3 

sat  4 

total 

pi 

p2 

p3 

p4 

total 

pi 

P2 

p3  1)4 

total 

run  1 

65 

87 

34 

70 

256 

0.15652 

0.1915 

0.1563 

0.1152 

0.1541 

5.3889 

8.083 

11.4  14.68 

9393 

run  2 

57 

74 

26 

56 

213 

0.30597 

0.3054 

0.2459 

0.2486 

03802 

5.4634 

8.529  937  1430 

9.424 

nin3 

59 

75 

24 

69 

227 

0.30435 

0.2910 

034 

03682 

03814  53414 

7.761 

9.83  11.6 

8.647 

run  4 

47 

75 

18 

51 

191 

0.30435 

03910 

0.24 

03682 

03814  5.3143 

6.945 

9.58  12.85 

8.667 

mnS 

67 

88 

25 

56 

236 

038 

0-1786  0.2059 

0.1579 

0.1993  5.8857 

9 

11.6  16-15 

10.14 

Mean 

59 

79.8 

25.4 

60.4 

224.6 

0.27024 

03515 

0.2176 

03116  03393  5.4587 

8.064 

10.3  13.90 

9394 

Var 

62 

49.7 

32.8 

73.3 

5963 

0,00416 

0.0037 

0.0014 

0.0050 

0.0035  0.0638 

0.609 

1.20  3.036 

0.409 

stdev 

7.8740 

7.0498 

5.7271 

8.5615 

24.4193 

0.06448 

0.0611 

0.0378 

0.0706 

0.0593  0  3526 

0.78 

1.10  1.742 

0.64 

Case3 

cum  chan  used  time  60 

#  good  trans 

#  StopG 

sat  1 

sat  2 

sat  3 

sat  4 

total 

pi 

p2 

p3 

p4 

total 

pi 

p2 

p3  p4 

total 

runl 

8 

8 

8 

8 

32 

104 

68 

5 

19 

196 

16 

13 

0  4 

33 

run  2 

8 

8 

8 

8 

32 

65 

86 

11 

14 

176 

8 

8 

6  5 

27 

run  3 

8 

8 

8 

8 

32 

76 

94 

13 

16 

199 

8 

16 

1  5 

30 

mn4 

8 

8 

8 

7 

31 

86 

76 

7 

21 

190 

10 

20 

2  1 

33 

mn5 

8 

8 

8 

8 

32 

89 

71 

7 

24 

191 

7 

16 

2  2 

27 

Mean 

8 

8 

8 

7.8 

313 

84 

79 

8.6 

18.8 

190.4 

9.8 

14.6 

2.2  3.4 

30 

Var 

0 

0 

0 

0.2 

03 

213.5 

117 

10.8 

15.7 

783 

133 

19.8 

5.2  3.3 

9 

stdev 

0 

0 

0 

0.4472 

0.44721 

14.6116 

10.817 

33863 

3.9623 

83487  3.6332  4.450  238  1.816 

3 
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Case3 

#stopNG 

gd/tot#msgs 

Time  to  transmit  gd  msgs 

sat  1 

sat  2 

sat  3 

sat  4 

total 

pi 

p2 

P3 

p4 

total 

pi 

p2 

p3  p4 

total 

run  1 

14 

75 

31 

85 

205 

0.7075 

0.3452  0.0877 

0.1027 

03345 

6.5673 

9.294 

9.8  14.16 

8332 

run  2 

15 

41 

34 

77 

167 

0.6373 

0.5 

0.1719 

0.08 

03431 

6.0154 

10.03 

13.5  13.93 

9.074 

run  3 

10 

50 

25 

94 

179 

0.7755 

0.4974 

0.2241 

0.0879 

03776 

6.0154 

10.03 

13.5  13.93 

9.074 

run  4 

7 

43 

36 

78 

164 

0.7611 

0.4222 

0.1207 

0.1235 

03647 

6.0154 

10.03 

13.5  13.93 

9.074 

runs 

12 

68 

27 

76 

183 

0.7607 

0.3817 

0.1273 

0.1341 

03557 

6.0154 

10.03 

13.5  13.93 

9.074 

Mean 

11.6 

55.4 

30.6 

82 

179.6 

0.7284 

0.4293 

0.1463 

0.1056 

03551 

6.1258 

9.887 

12.7  13,97 

8.925 

Var 

10.3 

233.3 

21.3 

57.5 

264.8 

0.0033 

0.0048 

0.0028 

0.0005 

0.0003 

0.0609 

0.11 

2.67  0.010 

0.11 

stdev 

3.2094 

15.274 

4.6152  7.5829 

16.2727 

0.0572 

0.069 

0.0528 

0.023 

0.0171 

0.2468 

0.331 

1.63  0.103 

0332 

Case4 

cum  chan  used  time  60 

#  good  trans 

#  stopG 

sat  1 

sat  2 

sat  3 

sat  4 

total 

pi 

P2 

P3 

p4 

total 

pi 

p2 

P3 

p4 

total 

run  1 

16 

16 

16 

16 

64 

81 

137 

33 

59 

310 

13 

12 

2 

11 

38 

run  2 

16 

16 

16 

16 

64 

94 

125 

35 

58 

312 

7 

22 

2 

11 

42 

run  3 

16 

16 

16 

16 

64 

89 

135 

34 

58 

316 

8 

13 

5 

8 

34 

run  4 

16 

15 

16 

16 

63 

97 

135 

37 

52 

321 

11 

6 

3 

10 

30 

runs 

16 

16 

16 

16 

64 

91 

129 

36* 

59 

315 

9 

23 

2 

8 

42 

Mean 

16 

15.8 

16 

16 

63.8 

90.4 

132.2 

35 

57.2 

314.8 

9.6 

15.2 

2.8 

9.6 

37.2 

Var 

0 

0.2 

0 

0 

0.2 

36.8 

25.2 

2.5 

8.7 

17.7 

5.8 

51.7 

1.7 

2.3 

27J 

stdev 

0 

0.4472 

0 

0 

0.44721 

6.0663 

5.0199 

1.5811 

2.9496 

4.2071 

2.4083 

7.190 

1.30 

1.517 

5.21536 

Ca$e4 

#stopNG 

gd/toti^m$gs 

Time  to  transmit  gd 

msgs 

sat  1 

sat  2 

sat  3 

sat  4 

total 

pi 

p2 

P3 

p4 

total 

pi 

p2 

P3 

p4 

total 

run  1 

14 

15 

2 

32 

63 

0.7043 

0.699 

0.6875 

0.3391 

0.5816 

5.8765 

8.036 

13.4 

14.41 

9.255 

run  2 

11 

13 

12 

50 

86 

0.7705 

0.6477 

0.5147  0.2959 

0.5389 

5.8765 

8.036 

13.4 

14.41 

9.255 

run  3 

11 

9 

6 

33 

59 

0.7417 

0.7337 

0.5763 

0.3295 

0.5863 

5.8765 

8.036 

13.4 

14.41 

9.255 

run  4 

11 

14 

2 

50 

77 

0.7405 

0.7542 

0.6271 

0.2889 

0.5847 

5.8765 

8.036 

13.4 

14.41 

9.255 

runs 

15 

10 

4 

43 

72 

0.6842 

0.6754 

0.6429 

0.3224 

0.5595 

5.8765 

8.036 

13.4 

14.41 

9.255 

Mean 

12.4 

12.2 

5.2 

41.6 

71.4 

0.7282 

0.702 

0.6097 

0.3152 

0.5702 

5.8765 

8,036 

13.4 

14.41 

9.255 

Var 

3.8 

6.7 

17.2 

77.3 

1173 

0.0012 

0.0018 

0.0044 

0.0005 

0.0004 

0 

0 

0 

0 

lE-14 

stdev 

1.9494 

2.5884 

4.1473 

8,792 

103305 

0.034 

0.043 

0.0663 

0.0217 

0.0206 

0 

0 

0 

0 

lE-07 
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APPENDIX  F.  ANALYSIS  OF  VARIANCE  RESULTS 


The  Analysis  of  Variance  results  for  the  model  are  shown  in  Tables  12-16. 


All  msgs 

#  good  trans 

#stopNG 

ANOVA 

ANOVA 

Source  of 

df 

SS 

MS  F 

Fcrit 

Source  of 

df 

SS 

MS  F 

Fcrit 

variation 

variation 

Jam 

1 

115368.1 

115368.1  265.9169 

4.494 

Jam 

1 

75645 

75645  178.0197 

4.494 

Chan 

1 

37584.45 

37584.45  S6.63006 

4.494 

Chan 

1 

30420 

30420  71,5891 

4.494 

JxC 

1 

7106.45 

7106.45  16.37997 

4.494 

JxC 

1 

4560.2 

4560.2  10.73178 

4.494 

Error 

16 

6941.6 

433.85 

Error 

16 

6798.8 

424.925 

Total 

19 

167000.6 

Total 

19 

117424 

#  stopG 

gd/tot#insgs 

ANOVA 

ANOVA 

Source  of 

df 

SS 

MS  F 

Fcrit 

Source  of 

df 

SS 

MS  F 

Fcrit 

variation 

variation 

Jam 

1 

2761.25 

2761.25  53.17766 

4.494 

Jam 

1 

0.369611 

0.369611  203.8603 

4.494 

Chan 

1 

312.05 

312.05  6.009629 

4.494 

•  Chan 

1 

0.121767 

0.121767  67.16105 

4.494 

JxC 

1 

2.45 

2.45i»CiiM 

4.494 

JxC 

1 

0.017422 

0.017422  9.609237 

4.494 

Error 

16 

830.8 

51.925 

Error 

16 

0.029009 

0.001813 

Total 

19 

3906.55 

Total 

19 

0.537809 

Time  to  transmit 

gd  msgs 

ANOVA 
Source  of 
variation 

df 

SS 

MS  F 

Fcrit 

Jam 

1 

0.030014 

0.030014S:I^g:: 

4.494 

Chan 

1 

0.64119 

0.641 1 9 

4.494 

JxC 

1 

0.004113 

0.0041 

4.494 

Error 

16 

2.489323 

0.155583 

Total 

19 

3.16464 

Table  18.  ANOVA  Results  for  All  Messages. 


Priority  1 

#  good  trans 

#stopNG 

ANOVA 

ANOVA 

Source  of 

df 

SS 

MS  F 

Fcrit 

Source  of 

df 

SS 

MS 

F 

Fcrit 

variation 

variation 

Jam 

1 

18361.8 

18361.8  193.5878 

4.494 

Jam 

1 

12954.05 

12954.05 

256.1354 

4.494 

Chan 

1 

336.2 

336.2  3.544544 

4.494 

Chan 

1 

61.25 

‘  61.25 

1.211073 

4.494 

JxC 

1 

16.2 

16.2  0.170796 

4.494 

JxC 

1 

92.45 

92.45 

L827978 

4.494 

Error 

16 

1517.6 

94.85 

Error 

16 

809.2 

50.575 

Total 

19 

20231.8 

Total 

19 

13916.95 

#  stopG 

gd/tot#insgs 

ANOVA 

ANOVA 

Source  of 

df 

SS 

MS  F 

Fcrit 

Source  of 

df 

SS 

MS 

F 

Fcrit 

variation 

variation 

Jam 

1 

288.8 

288.8  34-48358 

4.494 

Jam 

1 

0.61568 

0.61568 

147,6016 

4.494 

Chan 

1 

3,2 

3.2  0.38209 

4.494 

Chan 

1 

0.149826 

0,149826 

35-9189 

4.494 

JxC 

1 

5 

5S»;S»i5 

4.494 

JxC 

1 

0.049587 

0.049587 

U. 88777 

4.494 

Error 

16 

134 

8.375 

Error 

16 

0.06674 

0.004171 

Total 

19 

431 

Total 

19 

0.881832 

Time  to  transmit  gd  msgs 

ANOVA 

Source  of 

df 

SS 

MS  F 

Fcrit 

variation 

Jam 

1 

0.826772 

0.826772::|^|^$^3' 

4.494 

Chan 

I 

0.339027 

0.339027|5j^.^. 

4.494 

JxC 

1 

0.000624 

0.000624' 

4,494 

Error 

16 

1.007549 

0.062972 

Total 

19 

2.173972 

Table  19.  Priority  1  ANOVA  Results. 
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Priority  2 


fr  good  trans 

#stopNG 

ANOVA 

ANOVA 

Source  of 

df 

SS 

MS  F  Fcrit 

Source  of 

df 

SS 

MS 

F 

Fcrit 

variation 

variation 

Jam 

1 

22916.45 

22916.45  202.2635  4.494 

Jam 

1 

12400.2 

12400.2 

103.1202 

4.494 

Chan 

1 

5216.45 

5216.45  46.04104  4.494 

Chan 

1 

3225.8 

3225.8 

26.82578 

4.494 

JxC 

1 

2184.05 

2184.05  19.2767  4.494 

JxC 

1 

1584.2 

1584,2 

13.17422 

4.494 

Error 

16 

1812.8 

113.3 

Error 

16 

1924 

120.25 

Total 

19 

32129.75 

Total 

19 

19134.2 

#  stopG 

gd/tot#msgs 

ANOVA 

ANOVA 

Source  of 

df 

SS 

MS  F  Fcrit 

Source  of 

df 

SS 

MS 

F 

Fcrit 

variation 

• 

variation 

Jam 

1 

252.05 

252.05  7.839813  4.494 

Jam 

1 

0.61568 

0,61568 

147.6016 

4.494 

Chan 

i 

0.05 

0.05  0.001555  4.494 

Chan 

1 

0.149826 

0.149826 

35.9189 

4.494 

JxC 

1 

2.45 

2.45H1^K  4.494 

JxC 

1 

0.049587 

0,049587 

11.88777 

4.494 

Error 

16 

514.4 

32.15 

Error 

16 

0.06674 

0.004171 

Total 

19 

768.95 

Total 

19 

0.881832 

Time  to  transmit  gd  msgs 

ANOVA 

Source  of 

df 

SS 

MS  F  Fcrit 

variation 

Jam 

1 

1.136027 

1.136027fe^Sli  4.494 

Chan 

1 

3.375063 

3.375063pS?$S  4.494 

JxC 

1 

5.290529 

5.290529^W®;  4.494 

Error  16  3.092208  0.193263 


[Total _ 19  12.89383 _ 

Table  20.  Priority  2  ANOVA  Results. 


Priority  3 

#  good  trans 

#stopNG 

ANOVA 

ANOVA 

Source  of 

df 

SS 

MS  F 

Fcrit 

Source  of 

df 

SS 

MS  F 

Fcrit 

variation 

variation 

Jam 

1 

470.45 

470.45  54.07471 

4.494 

Jam 

1 

387.2 

387.2  14.4343 

4.494 

Chan 

1 

1110.05 

1110.05  127J92 

4.494 

Chan 

1 

980 

980  36.53308 

4.494 

JxC 

1 

661.25 

661.25  76.00575 

4.494 

JxC 

1 

649.8 

649.8  24.22367 

4.494 

Error 

16 

139.2 

8.7 

Error 

16 

429.2 

26.825 

Total 

19 

2380.95 

Total 

19 

2446.2 

#  stopG 

gd/tot#msgs 

ANOVA 

ANOVA 

Source  of 

df 

SS 

MS  F 

Fcrit 

Source  of 

df 

SS 

MS  F 

Fcrit 

variation 

variation 

Jam 

1 

76.05 

76.05  19.13208 

4.494 

Jam 

1 

0.165108 

0.165108  56.61683 

4.494 

Chan 

1 

11.25 

11-25  2.830189 

4.494 

Chan 

1 

0.32001 

0.32001  109.7342 

4.494 

JxC 

i 

4.05 

4-05^^|§|| 

4.494 

JxC 

1 

0.221268 

0.221268  75,87479 

4.494 

Error 

16 

63.6 

3.975 

Error 

16 

0.04666 

0.002916 

Total 

19 

154.95 

Total 

19 

0.753045 

Time  to  transmit  gd  msgs 

ANOVA 

Source  of 

df 

SS 

MS  F 

Fcrit 

variation 

Jam 

1 

9.40052 

9.40052^^^^ 

4.494 

Chan 

1 

5.072556 

5.072556^^^ 

4.494 

JxC 

1 

14.07058 

14.07058^^^; 

4.494 

Error 

16 

31.61032 

1.975645 

Total 

19 

60.15397 

Table  21.  Priority  3  ANOVA  Results. 
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Priority  4 
#  good  trans 

#stopNG 

ANOVA 

ANOVA 

Source  of 

df 

SS 

MS  F 

Fcrit 

Source  of 

df 

SS 

MS 

F 

Fcrit 

variation 

variation 

Jam 

1 

966.05 

966.05  IS,  17592 

4.494 

Jam 

1 

911.25 

911.25 

13.2401 

4.494 

Chan 

i 

4898.45 

4898.45  92.16275 

4.494 

Chan 

1 

6160.05 

6160.05 

89.50309 

4.494 

JxC 

1 

252.05 

252.05  4.742239 

4.494 

JxC 

1 

140.45 

140.45 

2.040683 

4.494 

Error 

16 

850.4 

53.15 

Error 

16 

1101.2 

68.825 

Total 

19 

6966.95 

Total 

19 

8312.95 

U  stopG 

gd/tot#nisgs 

ANOVA 

ANOVA 

Source  of 

df 

SS 

MS  F 

Fcrit 

Source  of 

df 

SS 

MS 

F 

Fcrit 

variation 

variation 

Jam 

1 

120.05 

120.05  10.48472 

4.494 

Jam 

1 

0.02551 

0.02551 

13.0608 

4.494 

Chan 

I 

162.45 

162.45  14.18777 

4.494 

Chan 

1 

0.157317 

0.157317 

80.54489 

4.494 

JxC 

1 

1.25 

4.494 

JxC 

1 

0.005166 

0.005166 

2.64S008 

4.494 

Error 

16 

183.2 

11.45 

Error 

16 

0.031251 

0.001953 

Total 

19 

466.95 

Total 

19 

0.219243 

Time  to  transmit  gdmsgs 

ANOVA 
Source  of 
variation 

df 

SS 

MS  F 

Fcrit 

Jam 

1 

0.096814 

0.096814Pt03^Ji 

4.494 

Chan 

1 

0.0195 

0.0195ppi^ 

4.494 

JxC 

1 

0.684107 

0.684107|l®iiffi 

4.494 

Error 

16 

53.80232 

3.362645 

Total 

19 

54.60274 

Table  22.  Priority  4  ANOVA  Results. 
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The  interaction  plots  for  the  number  of  messages  “stopG”  ANOVA  are  shown  in 


Figures  108-112. 


Figure  110.  Priority  2. 


Figure  112.  Priority  4. 


Figure  111.  Priority  3. 


5 


of  messages  #  of  messages 
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